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Nomenclature  of  Terms  Not  Defined  in  Text 


k  -_  Effective  multiplication  factor, 

ef  f 

e  Fast  fission  factor. 

p  Resonance  escape  probability. 

T]  Neutrons  born  per  thermal  absorption  in  fuel. 

f  Thermal  utilization. 

2  '         -2 

B  Buckling,  cm 

L  Thermal  diffusion  length,  cm. 

K  1/L,  cm 

2 
•£  Fermi  age,  cm  . 

v  Neutron  velocity,  cm/sec. 

q  Heat  generation  rate  in  center  of  fuel  at  start  of  pulse, 

°°  btu/hrft3. 

%  Reciprocal  period  of  power  rise,  sec 

A  ,X,         Parameters  for  empirical  fit  of  heat  generation  rate, 
btu/hrft3  and  sec-*  respectively. 

B  ,p  Parameters  for  empirical  fit  of  fuel  plate  surface  boundary 

condition,  °F  and  sec"*  respectively. 

L,L.  Thickness  of  fuel  and  moderator,  respectively  in  slab 

geometry,  cm. 

R,R.  Thickiess  of  fuel  and  moderator,  respectively  in  cylindrical 

geonv  try,  cm. 

x  Distance  in  fuel  from  the  cent  r  of  unit  cell,  slab  geometry, 

cm. 

x,  Distance  in  moderator  from  outside  of  unit  cell,  slab 

geometry,  cm. 

r  Distance  An   fuel  and  moderator  from  center  of  unit  cell, 

cylindrical  geometry,  cm. 

i,ji,n         Summation  indices. 


f,m  Subscripts  denoting  fuel  and  moderator,  respectively. 

p,s  No.  of  terms  in  heat  generation  rate  and  interface 

boundary  condition  approximations. 


1  .0  INTRODUCTION 

The  concept  of  reactor  safety  Is  extremely  important  in  the  engineering 
application  of  nuclear  power  systems.   The  United  States  Atomic  Energy 
Commission  has  therefore  authorized  an  extensive  study  in  this  area. 
This  investigation  uses  experimental  data  resulting  from  that  study  to 
attempt  to  define  the  mechanisms  of  reactor  shutdown. 

The  safety  of  a  nuclear  reactor  system  is  usually  considered  in  terms 
of  its  void  and  temperature  coefficients  of  reactivity.   If  one   designs 
a  reactor  in  such  a  manner  as  to  make  both  of  these  coefficients  negative, 
the  system  will  tend  to  stabilize  itself  if  some  external  perturbation 
is  placed  on  the  system.   This  is  due  to  the  fact  that  when  the  excess 
reactivity,  k  ff-l,    of  a  system  is  increased,  the  power  tends  to  rise, 
thereby  increasing  the  temperature  of  the  system,  and  this  in  turn 
causes  a  decrease  in  reactivity.   In  the  case  of  a  liquid  moderated 
system,  voids  may  be  introduced  whic  i  will  further  decrease  the  reactivity. 

If  one  considers  the  bare  reactor  age-diffusion  criticalitv  equation, 

k    =  £pnf  e 
eff — -z ,  it  is  possible  to  see  how  these  effects  manifest 

1  +  L  B 


themselves  in  the  nuclear  constants  (22).   In  the  case  of  the  void 
coefficient,  the  significant  effect  is  on  the  Fermi  age,  j£.   Since  part 
of  the  moderator  is  removed,  the  age  increases,  thus  decreasing  the  fast 

non-leakage  factor,  e .   In  the  case  of  the  temperature  coefficient 

there  are  several  tf feces  thac  must  be  considered.   First,  the  moderator 
expands  because  of  the  positive  coefficient  of  expansion  of  the  liquid. 


This  results  in  a  decrease  in  the  density  of  the  moderator,  and  if  the 
reactor  is  under-moderated  the  fast  non-leakage  factor  and  p  will  decrease. 
Thus  it  is  evident  that  one  safety  criteria  is  that  the  core  should  always 
be  slightly  under-moderated  even  though  this  will  increase  the  critical 
mass.   Second,  the  fuel  elements  expand,  expelling  more  of  the  moderator 
from  the  core  causing  k  ff  to  decrease;  however,  simultaneously,  the 

effective  size  of  the  core  increases,  causing  a  decrease  in  the  buckling, 

2 
B  ,  of  the  system.   This  decrease  in  buckling  has  both  a  positive  and 

a  negative  effects  upon  reactivity  or  k  rc.      The  increase  in  size  increases 

eff 

-B  v 
the  non-leakage  factor  for  both  fast  and  thermal  leakage,  e and 

2  2 
1/(1  +  L  B  ),  respectively.   The  removal  of  moderator  increases  ^ 

reducing  reactivity  as  described  in  the  void  coefficient  discussion.   All 

of  the  effects  described  above   with  the  exceptioi  of  the  buckling, 

must  be  considered  to  be  of  a  delayed  nature. 

Another  group- of  effects  exists  which  affect  the  reactivity 

immediately  and  these  effects  are  therefore  classified  as  prompt  (41). 

The  first  of  these  prompt  effects  is  Doppler  broadening  (38).   Because 

of  the  increased  kinetic  energy  of  U    target  nuclei  with  increased 

temperature,  the  width  of  resonance  absorption  is  increased,  but  the 

height  of  the  peak  is  decreased  (22),  because  the  total  area  beneath 

the  resonance  curve  remains  constant   If  the  resonance  absorption  cross 

sections  are  large,  so  that  essentially  all  neutrons  with  energies  in  the 

resonance  region  are  captured,  the  widening  of  the  region  will  result  in 

a  decrease  in  the  resonance  escape  probability,  p,  and  thus  the  reactivity 

decreases  as  the  temperature  rises.   A  second  prompt  effect  is  caused  by 

hardening  of  the  thermal  neutron  spectrum  as  the  moderator  temperature 


rises.   The  result  of  this  hardened  spectrum  is  that  the  average  thermal- 

2 
neutron  velocity  increases,  _L_  increases  and  the  thermal  non-leakage 

1 
factor,  — -   ,  decreases.   One  would  expect  that  the  neutrons  per 

1  +  LB 
absorption  would  also  be  affected  by  this  spectral  change  but  assuming 

the  normal  1/v  dependence  for  all  cross  sections,  this  effect  cancels 

since  the  terms  comprising  this  effect  are  composed  of  the  ratios  of 

cross  sections.  '  One  final  effect  which  is  neither  part  of  the  temperature 

coefficient  nor  the  void  coefficient  must  be  considered.   This  is  the 

formation  of  radiolytic  gases. 

The  reasons  that  one  is  inclined  to  speak  of  prompt  and  delayed 
coefficients  of  reactivity  is  that  for  power  bursts  of  low  reactivity 
and  correspondingly  long  periods,  the  delayed  effects  may  play  a  great 
part  in  the  shut-down  mechanism.   However,  if  one  supposes  a  burst  with 
a  very  short  period  then  it  is  obvious  that  these  delayed  effects  will 
not  have  had  time'  to  act.   What  constitutes  a  "short  time"  can  be  answered 
by  determining  the  heat  transfer  time  constant  of  the  fuel  elements. 
This,  in  part,  is  the  subject  of  the  proposed  investigation.   Iriarte  (27) 
reported  heat  transfer  time  constants  for  infinitely  long  cylindrical 
UO   fuel  elements  surrounded  by  a  helium  film  which  sei'ved  as  a  thermal 
1  ond  and  which  were  clad  with  zirconium,  stainless  steel  or  aluminum. 
While  the  data  presented  by  Iriarte  are  not  applicable  to  the  SPERT-I  or 
TRIGA  Reactors,  the  techniques  may  be  useful  in  determining  the  relative 
importance  of  the  delayed  effects. 

Before  considering  further  the  scope  of  this  project  it  is  informative 
to  investigate  the  prior  work  in  the  field.   During  the  early  summer  of 
1954  a  series  of  experiments  were  made  on  the  BORAX-I  Reactor  to  investigate 
the  ability  of  the  reactor,  when  opeiated  in  the  subcooled  (non-boiling) 


condition,  to  protect  itself  against  the  results  of  sudden,  artificially 

induced  increases  in  reactivity.   Inasmuch  as  this  set  of  experiments 

completed  the  program  for  the  BORAX-I  Reactor,  the  final  runaway  experiment 

was  Intentionally  made  under  conditions  which  led  to  destruction  of  the 

reactor.   In  the  final  experiment,  a  control  rod  worth  four  per  cent 

k  __,  was  ejected  from  the  reactor  core,  inducing  an  exponential  power 
erf 

increase  which  had  a  period  of  2.6  milliseconds.   This  final  experiment 
resulted  in  a  melting  of  most  of  the  fuel  plates  and  failure  of  the 
reactor  tank.   Fuel  plate  fragments  were  scattered  for  a  distance  of  200 
to  300  feet  (8).   This  set  of  experiments,  along  with  the  earlier  operation 
of  the  BORAX-I,  (7),  established  two  important  safety  axioms  for  water 
moderated  reactors.   First,  for  any  given  system  there  is  a  reactivity 
insertion  beyond  which  the  reactor  cannot  react  fast  enough  to  shut  itself 
down  before  damage  is  done,  and  second,  water  moderated  reactors  can  be 
designed  to  have  a  high  degree  of  inherent  self -protection  against  the 
effects  of  sudden  large  reactivity  increases.   A  less  critical  fact  that 
resulted  from  these  tests  was  that  if  the  transients  are  started  at 
boiling  conditions  (such  as  in  a  boiling  water  reactor),  the  maximum 
power  and  fuel-plate  temperature  reached  are  less  than  if  the  transient 
is  started  with  the  reactor  in  a  subcooled  condition.   This  is  as  would 
be  expected  if  void  formation  due  to  boiling  were  the  shut-down  mechanism, 
since  for  the  subcooled  system  the  reactor  could  actually  achieve  a 
stable  positive  period  before  any  negative  reactivity  would  result. 
This,  along  with  the  observation  of  large  quantities  of  steam  and  water 
being  expelled  from  the  system.,  led  to  the  conclusion  that  it  was  the 
void  formation  which  was  shutting  the  reactor  down.   Vvith  this  background 
in  mind,  the  Atomic  Energy  Commission  set  out  on  an  intensified  program 


to  determine  empirically  the  safe  upper  operating  limits  on  each  class 
of  reactors  including  the  pressurized  and  unpressurized ,  boiling  and 
non-boiling  thermal  reactors,  both  heterogeneous  and  homogeneous,  and 
fast-reactor  systems  (23).   This  program  initiated  the  SPERT  (Special 
Power  Excursion  Reactor  Tests)  and  KEWB  (Kinetic  Experiments  on  Water 
Boilers)  programs.   SPERT  is  the  heterogeneous  reactor  test  facility 
and  KEWB  is  the  design  for  the  same  type  of  tests  on  a  homogeneous  reactor 
test  facility.   In  1956,  W.  B.  Nyer,  et .  al.  (39),  reported  the  results 
of  the  initial  transient  test  on  the  SPERT-I  facility  and  concluded 
that  the  SPERT-I  (43)  reactor  demonstrated  qualitatively  the  results  of 
the  BORAX-I  experiments  although  SPFRT-I  was  more  stable  after  the 
initial  power  burst.   Factors  which  could  contribute  to  this  difference 
In  behavior  were  the  known  differences  in  the  fuel  assembly  construction, 
the  possible  differences  in  the  effective  void  coefficients,  and/or 
the  differences  in  the  reflector-tank  environment.   Another  important 
point  discovered  at  this  time  was  that  there  was  a  unique  relationship 
between  the  peak  power  and  the  transient  period.   The  data  were  fitted 
rather  well  by  two  straight  lines  on  a  plot  of  log  reciprocal  period 
versus  peak  power.   The  slopes  of  these  line  were  approximately  0.8  and 
1.7  for  the  lower  and  upper  regions  respectively.   The  point  of  inter- 
section  occurred  at  Ak  =  0.74%,  or  ab  jut  prompt  critical.   In  June  of 
1957,  R.  W.  Miller  reported  on  some  interesting  work  in  analyzing  the 
reactivity  behavior  during  SPERT-I  transients  (34).   Miller  pointed  out 
that  it  was  not  necessary  to  remove  all  of  the  initially  inserted  reactivity 
to  limit  a  power  excursion.   This  is  a  result  of  the  delayed  neutrons, 
which  for  a  very  short  period  excursion,  do  not  contribute  to  the  flux 
during  the  rise  in  power.   Thus,  for  short  excursions  the  reactivity 


compensation  at  maximum  power  need  only  he  Ak(l-fO-P.   Ak  is  the  initial 
change  in  k  ff  and  (3  is  the  delayed  neutron  fraction.   The  total  Ak  must 
be  accounted  for  in  very  long  period  transients.   In  the  case  of  inter- 
mediate period  transients  the  compensated  reactivity  was  calculated  as 
a  function  of  the  reciprocal  period  by  numerical  solution  of  the  reactor 
kinetics  equations  using  the  experimentally  determined  power  traces. 
Later  in  1957,  G.  0.  Bright,  et.  al.,  suggested  a  model  for  reactor 
burst  behavior  (2),  based  on  the  earlier  work  of  Klaus  Fuchs  at  Los 
Alamos  (21).   This  model  postulated  a  shutdown  effect  proportional  to 
the  energy  release.   However,  the  model  provided  no  method  for  the  shu*  - 
down  energy  to  be  removed  and  allowed  for  no  time-delays  between  the 
energy  release  and  the  appearance  of  the  shutdown  effect.   Much  later 
this  model  was  further  modified  by  S.  G.  Forbes  (13)  who  let  the  shut- 
down effect  be  proportional  to  the  energy  release  raised  to  some  power, 
n,  and  allowed  for  some  arbitrary  delay  time.   In  this  analysis  the 
model  was  tested  against  some  experimental  transients  from  SPERT-I  and 
values  of  n  from  1.5  to  2  were  successively  used.   The  results  also 
showed  that  the  delay  time  was  significant  in  matching  the  data;  however 
no  exact  value  of  the  delay  time  was  determined.   The  overall  effect  of 
these  works  is  to  convince  one  that  the  primary  shutdown  mechanisms 
are  intimately  tied  up  with  the  energy  release  although  no  real  infor- 
mation on  the  exact  phenomena  can  be  found.   In  1958,  Griff ing  and 
Deverall  (10)  coupled  the  energy  shut-down  model  with  the  reactor  kinetics 
equation  including  six  delay  groups  and  again  showed  that  the  energy  model 
could  describe  qualitatively  the  power  traces  obtained  in.  the  SPERT-I 
transients  even  long  after  the  initial  burst.   This  work  used  a  mathe- 
matical structure  of  the  shut-down  equation  considered  much  earlier  (1951) 


by  Chernick  (4)  for  no  delay  groups  and  in  1956  by  Margulles  (31)  with 
one  delay  group. 

In  December  of  1958,  Deverall  and  Griff ing  (9)  reported  the  first 
attempts  at  trying  to  relate  the  shutdown  reactivity  with  the  thermo- 
dynamic characteristics  of  the  SPERT-I  system.   In  this  report  the  change 
in  reactivity  for  transients  with  long  periods  was  related  to  the  temperature 
rise  in  the  moderator  alone  although  the  fact  was  recognized  that  the 
fuel  element  temperature  -  ise  should  also  be  considered.   Considering 
only  the  moderator  temperature  rise,  they  found  reactivity  compensations 
at  peak  power  that  were  roughly  one-half  of  those  reported  by  Miller  (34). 
Since  they  felt  that  the  data  with  which  they  were  working  were  only 
accurate  to  within  a  factor  of  t--o,  they  did  not  pursue  the  investigation 
further.   In  January  of  1958,  Horning  of  Ramo-Woolridge  reported  on  a  model 
for  transients  in  SPERT-I  (18).   This  report  develops  a  general  model, 
taking  into  account  the  void  formation  and  the  thermal  expansion.  However, 
no  real  attempt  was  made  to  interpret  these  constants  in  temis  of  the 
distribution  of  energy  in  the  fuel  and  moderator  or  he  nuclear  and  thermo- 
dyn-'mJc  const  nts  of  the  system.    nother  report, bo  nd  ander  the  same 
cover,  by  H.  C.  Corben  (26)  tirats  the  problem  of  oscillations  found  after 
the  burst  as  the  power  approaches  some  stead}  s  ate  level.   Although  the 
mechanism  responsible  for  the  oscillations  need  not  be  the  same  as  the 
shut-down  mechanism,  there  is  certainly  the  possibility  that  they  are  one 
and  the  same.   A  third  report  by  3.  Bi: khof i  (26)  treats  the  problem  of 
void  formation  from  the  r ji  '  of  view  of  the  growth  of  bubbles.   The 
"Bubble  Void"  is  probaDly  the  dominant  shutdown  effect  in  a  certain  type 
of  excursion  such  as  in  the  initial  BORAX-I  experiments.   The  analytic 
representation  of  this  effect  is  the  least  known  and  thus  is  being  widely 


sought.   An  extensive  study  of  bubble  formation  including  a  critical 
review  of  the  literature,  an  evaluation  of  the  merits  of  purely 
theoretical  approaches  to  the  development  of  a  void  model,  and  an 
investigation  of  the  possible  formulation  of  the  nucleate  boiling  void 
was  reported  by  the  Vitro  Engineering  Company  in  May  of  1959  (28). 

Although  an  exact  definition  of  bubble  formation  may  not  be  within 
the  scope  of  thi,s  investigation,  the  determination  of  the  heat  flow  into 
the  moderator  and  the  transient  temperature  distribution  in  the  moderator 
should  shed  some  light  on  even  this  difficult  problem. 

In  July  of  1958,  J.  C.  Haire  (25)  reported  the  results  of  a  great 
number  of  the  SPERT-I  transients.   This  report  presented  data  on  the 
reactor  power,  fuel  plate  surface  temperatures  and  pressures  as  a  function 
of  time  during  the  transients.   These  are  the  data  that  will  be  used 
extensively  in  the  initial  phases  of  the  investigation  proposed  herein. 

During  late  1958  and  1959,  several  models  were  proposed  to  investigate 
and  explain  the  inherent  shutdown  characteristics  of  the  SPERT-I  reactor. 
The  "Conduction  Boiling"  model  suggested  by  S.  G.  Forbes  (14)  is  certainly 
credible  in  that  it  takes  into  account  the  flow  of  heat  into  the  moderator 
in  a  much  more  exact  manner  than  any  of  the  earlier  investigations.   This 
model  was  quite  successful  in  predicting  the  power,  energy  release  and 
temperature  at  the  time  of  peak  power  as  a  function  of  the  reciprocal 
Deriod.   However  the  model  still  represents  the  shutdown  mechanism  in  terms 
of  empirical  parameters.   Also  the  non-boiling  shutdown  effects  are  not 
taken  into  account.   The  "Clipped  Exponential"  model  suggested  by  R.  W. 
Miller  (35)  made  some  very  useful  assumptions  on  the  shape  of  the  reactor 
Dower  burst  to  ease  the  analytical  solution  of  kinetics  equations.   While 
this  model  produced  some  useful  criteria  in  the  understanding  of  self-shutdown 


it  again  made  use  of  lumped  parameters  which  were  not  easily  interpreted 
in  terms  of  the  thermodynamic  and  nuclear  characteristics  of  the  system. 
E.  T.  Clark  (5,6)  as  early  as  1956  had  postulated  a  prompt  fission  pro- 
duct having  a  large  absorption  cross  section  for  thermal  neutrons  as  an 
explanation  of  the  self-shutdown  of  power  excursions  observed  in  the 
SPERT-I  Reactor.   Later  evidence  (42)  seems  to  indicate  that  this  model 
is  less  likely  to  be  valid  than  the  more  conventional  models. 

Also  in  1958,  General  Atomic  introduced  their  TRIGA  Reactor  with 
zirconium  hydride  moderator  which  demonstrated  a  larger  prompt  shutdown 
mechanism  than  either  the  BORAX-I  or  the  SPERT-I  (41).   In  this  case  the 
reactor  was  designed  to  have  a  shutdown  mechanism  which  would  act  by 
hardening  the  thermal  neutron  energy  spectrum  thus  increasing  the  thermal 
leakage  to  cause  shutdown.   The  reason  for  choosing  this  effect  was  that 
they  felt  that  it  would  act  more  quickly  and  thus  provide  a  safer  reactor 
than  the  accepted' moderator  expansion  and  expulsion  mechanisms.   The 
spectrum  effect,  so  important  in  the  TRIGA  Reactors  must  also  act  to  some 
extent  in  SPERT-I.   The  amount  of  this  effect  has  apparently  never  been 
determined . 

P.  French,  in  1959,  (18)  reported  an  attempted  solution  of  the  transient 
heat  conduction  equation  in  the  fuel  and  moderator  to  determine  the 
temperature  distribution.   However,  this  work  appears  to  be  in  error  in 
that  he  forced  a  separation  of  variables  solution  on  the  equations  whereas 
the  spatial  and  time  dependence  cannot  be  expressed  as  a  simple  product 
except  after  sufficiently  long  times  so  that  the  transient  term  has  dis- 
appeared.  H.  L.  McMurry  (32)  reported  the  temperature  distribution  in  a 
fuel  plate,  cladding  and  moderator  with  exponentially  rising  power  for 
pure  conduction.   This  is  an  excellent  piece  of  work  but  the  mathematical 
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model  turns  out  to  be  more  difficult  than  is  either  warranted  or 
necessary  for  the  analysis  of  the  SPERT-I  data.   The  McMurry  report  makes 
mention  of  earlier  work  by  H.  Greenspan  (24)  on  the  same  problem  with 
similar  results.   Several  investigators,  Kattwinkel  (29),  Kirchenmayer 
(30),  Stein  (40),  Epel  (11),  Arpaci  and  Clark  (1),  Ermakov  and  Ivanov  (12), 
report  analytical  solutions  to  the  transient  heat  conduction  equation. 
However,  none  of ( them  were  working  on  the  problem  with  reference  to  the 
SPERT-I  investigations  and  their  results  are  not  directly  applicable  to 
the  use  of  the  SPERT-I  experimental  measurements.   As  a  result  the 
temperature  distribution  in  the  fuel  and  moderator  during  a  transient  is 
not  well  known  to  date. 

In  July  of  1959,  Forbes,  et.  al.  (16),  summarized  the  work  done  up 
to  that  time.   They  showed  that  by  fitting  empirically  the  "Conduction 
Boiling"  model  to  the  SPERT-I  data  and  including  the  effect  of  moderator 
and  fuel  element  expansion  they  could  fit  the  experimental  compensated 
reactivity  at  peak  pow^r  versus  reciprocal  period  curves  for  values  of 
the  reciprocal  period  greater  than  5  seconds   .   They  also  showed  that 
for  values  of  the  reciprocal  period  greater  than  20  seconds   ,  the  steam 
void  contribution  to  reactivity  was  considerable.   For  the  longer  period 
region  they  postulated  an  additional  shutdown  effect  from  radiolytic  gases. 

This,  of  course,  lends  considerable  credence  to  the  "Boiling  Conduction" 
model,  however  it  does  have  one  glaring  shortcoming.   To  extend  it  to 
another  reactor  requires  at  least  one  transient  burst  experiment  to 
determine  the  parameters.   However,  if  the  fundamental  mechanism  were 
understood  exclusively  in  terms  of  the  basic  nuclear,  thermodynamic  and 
hydrodynamic  characteristics  of  the  system,  one  could  confidently  predict 
the  limits  of  safe  ope re:  Lona  for  different  reactor  systems. 
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The  satisfactory  application  of  the  "Boiling  Conduction"  model 
led  to  two  sets  of  experiments  designed  to  test  that  model  and  the 
postulated  radiolytic  gas  effect.   The  first  set  of  experiments  (20) 
showed  with  reasonable  certainty  that  radiolytic  gas  formation  was  not 
a  primary  contributor  to  self-shutdown  in  the  SPERT-I  Reactor.   The 
second  set  of  experiments  (19)  consisted  of  coating  all  of  the  fuel  plates 
with  approximately  five  mils  of  insulation,  Lithocote  LC-34,  and  of 
performing  transient  tests  on  the  reactor.   Some  of  the  tests  were  run 
with  transients  of  such  magnitude  that  boiling  would  occur  in  the  bare 
core  and  not  in  the  insulated  core.   The  remaining  tests  involved 
transients  in  which  there  would  be  boiling  in  both  cores,  and  the 
differences  in  heat  transfer  rates  were  expected  to  be  reflected  as 
changes  in  the  reactor  behavior.   Power  burst  shapes  for  transients  of 
the  same  period  in  both  the  bare  and  insulated  cores  were  essentially 
the  same.   In  view  of  the  identical  reactor  behavior  for  the  bare  and 
insulated  tests,  it  would  appear  that  the  core  insulation  produced  no 
appreciable  effects  on  the  shutdown  mechanism.   Since  it  seems  reasonable 
to  assume  that  any  shutdown  effect  due  to  boil  ins  would  be  effected  by 
the  core  insulation,  boiling  would  not  seem  to  play  any  part  in  the 
self-shutdown  mechanism.   However,  if  the  heat  transfer  rate  was  small 
enough  there  would  be  negligible  temperature  drop  across  the  insulation 
and  thus  the  effect  on  boiling  might  be  unchanged  by  the  insulation.  This 
experiment  shows  clearly  the  need  for  detailed  calculation  of  the  temperature 
distributions  in  the  fuel  and  moderator  during  the  transients.   One  final 
report  should  be  cited  in  this  summary.   In  April  of  1960,  Miller  (36) 
reported  on  some  photographic  investigation^  of  boiling  during  transients 
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in  SPERT-I.   These  results  clearly  Indicated  that  boiling, was  an  important 
agent  in  the  initial  reactor  self -shutdown  whenever  the  fuel  plate 
temperature  was  sufficiently  high. 


13 


2.0  THEORY 
2.1  Derivation  of  Equations 

The  direct  approach  to  determining  the  s e If -shutdown  effects  must 
include  the  determination  of  the  temoerature  distributions  in  the  fuel 
elements  and  moderator  throughout  the  core  during  a  transient.   This 
problem  can.  be  accomplished  by  investigating  the  exact  temperature 
distribution  in  a  center  element  only  and  relating  all  effects  to  this 
center  element.   The  exact  solution  of  the  multi-region,  transient 
heat  conduction  equations  even  for  a  single  fuel  element,  making  use  of 
only  the  power  versus  time  data  from  the  SPERT-I  transients,  is 
exceedingly  difficult  as  pointed  out  in  the  work  by  McMurry  (32).   How- 
ever, if  use  is  made  of  the  available  fuel  plate  surface  temperature  data 
as  well  as  the  oower  data  during  a  transient,  the  problem  is  reduced  to 
two  single  region  problems.   Although  the  solutions  are  much  simpli- 
fied over  the  two  region  problem,  they  are  still  complicated  and  therefore 
have  been  programmed  for  the  Kansas  State  University  IBM  650  computer. 

The  methods  for  determining  the  steady  state  temperature  distribution 
throughout  a  unit  cell  of  thermal,  heterogeneous  liquid-moderated  reactor 
are  discussed  in  Nuclear  Engineering  (41)  by  C .  F.  Bonilla  and  in 
Nuclear  Reactor  Physics  (37)  by  R.  L.  Murray.   It  will  be  considered 
sufficient  for  this  work  to  outline  the  differences  that  must  be  accounted 
for  in  transient  operation. 

The  partial  differential  equation  for  conductive  heat  transfer 
applicable  in  the  fuel  and  moderator  of  a  nuclear  power  reactor  during 
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a  transient  but  before  boiling  is  established  is 

V26(x,t)  +S*5it)   =  i/a-S9i£ttl  (1) 

k  3t 

2 
where  V   is  the  Laplacian  operator  (44),  9  is  the  temperature  rise  above 

the  initial  temperature  [T(x,t)-T  ],  x  is  the  position  variable,  k  is 
the  mean  thermal  conductivity,  a  is  the  mean  thermal  diffusivity,  t  is 
the  time  variable  and  £  is  the  volumetric  heat  generation  rate.   This 
equation  in  slab  and  cylindrical  geometries  is  directly  applicable  to 
the  analysis  of  transient  behavior  in  a  reactor  following  a  change  in 
reactivity.   The  derivation  of  solutions  to  this  equation  in  the  fuel 
and  moderator  during  power  transients  will  be  shown  in  detail  for  slab 
geometry  (Appendix  A)  and  the  important  elements  of  the  solution  in 
cylindrical  geometry  will  be  tabulated  in  Section  2.2.1. 

In  any  nuclear  reactor  there  are  heating  effects  due  directly  to 
fission  fragments  and  to  the  attenuation  of  other  nuclear  particles. 
Within  the  fuel  element,  fission  heating  far  overrides  the  other  at- 
tenuation effects.   Therefore,  the  heat  generation  rate,  qf(x,t),  in  the 
fuel  elements  is  proportional  to  the  thermal  neutron  flux.   The  neutron 
flux  during  a  transient  can  be  expressed  as  a  simple  product  of  the 
spatial  and  the  time  dependencies.   Thus  the  heat  generation  rate  is 
also  separable  in  space  and  time  as  shown  in  equation  (2). 

qf(x,t)  =  ff(x)   g  (t)  (2) 

The  spatial  dependence,  ff(x),  can  be  obtained  easily  from  the  steady 
state  analysis  (40)  and  is  given  in  equation  (3). 

f£(x)  =  q   cosh  kx  (3) 

t       oo 
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Here  q    is  the  heat  generation  rate  at  the  center  of  the  fuel,  x   is 

oo  .  — 

the  inverse  thermal  neutron  diffusion  length  and  x  is  the  distance  from 

the  center  of  the  fuel. 

The  time  dependence,  gf(t),  of  the  neutron  flux  and  thus  the  heat 

eneration  rate  is  expressed  as  the  sum  of  exponentials,  thus 

s 
gf(t)  =  £  a.eY  (A) 

Substituting  equations  (3)  and  (4)  into  equation  (2),  an  expression  for 


qf(x,t)  is  obtained;  that  is 


qf(x,t)  =£   q00cosh(Kx)a.e  j 
J=l 


^  (5) 


Substituting  equation  (5)  into  equation  (1)  yields  the  differential 
equations  which  must  be  solved  to  obtain  the  temperature  distribution 

in  the  fuel,  that  is 

s 

r^2D  i      *.\  l  V  q  cosh(Kx)a.e  j    _  ,.   aGL(x,t)  ,,» 

V  Or(x,t)  +  )   Moo i   J    =  1/a   f  (6) 

,_1      k  at 

Now  from  investigation  of  the  heat  generation  rate,  q  (k,t)  in  the 

moderator,  it  is  noted  that  there  is  no  heating  due  to  fission,  but  there 

is  heating  due  to  nuclear  particles  which  stream  out  of  the  fuel  and  are 

attenuated  in  the  moderator.   The  moderator  heating  is  approximately  5 

to  1%   (3)  of  the  recoverable  energy  from  fission  and  is  sufficiently 

uniform  in  space  to  be  so  considered.   Therefore,  q  (x,t)  in  the  moderator 

is  indepenuent  of  the  spatial  variable  but  still  is  time  dependent  as 

shown  in  quation  (7). 

s 

q  (x,t)  =V   F  a.eV  C7) 
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The  heating  effects  in  the  moderator  are  proportional  to  the  neutron 
flux  in  the  fuel;  thus  F  is  the  fraction  of  the  recoverable  energy  released 
in  the  fuel  which  is  dissipated  in  the  moderator. 

Substitution  of  equation  (7)  into  equation  (1)  yields  the  differen- 
tial equation  to  be  solved  for  the  temperature  distribution  in  the  moderator, 

that  is 

s         x   t 

r-72a  (      -%   V   F  a,e  j   _  ,,   a  0  (x,t)  ,Qv 

V  0  (x,t)  +  >      i   J   =  1/a    m  (8) 

r,       k  at 

J=.l 

Equations  (6)  and  (8)  form  the  general  set  which  must  be  solved 
to  obtain  the  temperature  distribution  in  a  unit  cell  of  the  reactor  during 
a  transient.   The  geometry  of  the  unit  cell  in  which  these  equations  must 
be  applied  is  shown  in  Figure  1.   The  simultaneous  solution  of  these 
equations  is  evtremely  difficult.   The  availability  of  the  fuel  element 
surface  temperature  as  a  function  of  time  during  the  transients  greatly 
simplifies  this  situation.   The  problem  is  reduced  to  solving  the  equations 
independently  in  the  fuel  and  moderator,  using  the  experimentally  measured 
temperatures  at  the  interface  as  a  boundary  condition  for  both  equations. 
The  other  boundary  condition  necessary  in  each  case  is  a  zero  heat  flow 
condition  at  the  center  of  the  unit  cell  for  the  fuel  regions  and  at  the 
boundary  of  the  unit  cell  for  the  moderator. 

2.2  Analytical  Solutions 

The  -time  dependent  thermal  diffusion  equations  in  the  fuel  and 
moderator  can  be  solved  for  the  temperature  distribution  assuming  that 
conduction  is  the  primary  mode  of  heat  transfer.   The  equation  will  hold 
for  all  time  in  the  fuel  plate.   In  an  attempt  to  represent  as  well  as 
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Figure  I.  Geometry  and  boundary   conditions  in  the  unit  cell  used 
to  determine  the  temperature  distributions. 
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possible  the  experimental  data  which  is  used  as  input  to  this  problem 
and  to  allow  some  flexibility  in  the  application  of  these  equations,  the 
differential  equations  were  solved  subject  to  several  forms  of  represen- 
tation of  the  boundary  conditions  and  the  forcing  functions.   In  one 
case  the  problem  was  solved  in  cylindrical  geometry. 

In  slab  geometry  the  one-dimensional  solutions  to  the  transient  heat 
transfer  equation,  in  which  the  time  dependence  of  heat  generation  rate 

and  the  fuel  element  surface  temperature  are  represented  empirically 
s  p 

ZX  t      V-1    8  t 
A  e  j   and  )   B  e  i  respectively,  are  derived  in  Appendix  A 

j=l  i=l 
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The  equivalent  solutions  in  cylindrical  geometry  are  derived  in 
Appendix  A  and  are  given  here  as 
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where  'J  's  are  tTie  roots  of  the  equation,  J  (x)  =  0 
n  o 

and  p  '  s  are  the  roots  of  the  equation, 
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The  solutions  to  the  transient  heat  transfer  equations  in  slab 

geometry  in  which  the  time  dependence  of  the  heat  generation  rate  and 

the  first  derivative  of  the  surface  temperature  with  respect  to  x  are 
s'  p* 
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in  the  fuel  and  moderator,  respectively. 

The  solution  to  the  transient  heat  transfer  equation  in  slab 

geometry  in  which  the  time  dependence  of  the  heat  generation  rate  and 

s  p 

the   surface   temperature  are  represented  by     )        A     e   j  and      )        B     cos  p   t 
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respectively  are  derived   in  Appendix  A  and   are  given  here  as 
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2.2.1  Temperature  Distributions.   Equations  (15)  and  (16)  of  the 

previous  section  can  be  evaluated  to  obtain  the  temperature  as  a  function 

of  position  and  time  in  any  unit  cell  of  a  reactor  if  the  heat  generation 

rate,  q(x,t),  and  the  fuel  surface  temperature,  9f(L,t),  are  known  and 

can  be  expressed  in  the  appropriate  analytical  form.   The  experimental 

values  of  these  variables  during  applicable  transient  tests  on  the 

SPERT-I  reactor  system  were  obtained  in  graphical  form  from  "Sub-cooled 

Transient  Tests  in  the  SPERT-I-A  Reactor  -  Experimental  Data"  by 

J.  C.  Haire  (25).   Numberical  power  and  temperature  data  were  obtained 

from  the  graphs.   These  numerical  data  were  normalized  to  a  zero  initial 

P 

temperature  then  fit  empirically  by  an  even  trigonometric  series,  )      A  cos  p.t, 

i=l 

for  the  temperature  traces.   The  power  traces  were  reduced  to  give  the  heat 

generation  rate  in  the  center  of  a  central  fuel  element  and  moderator 

s 

V     X.  t 
region  and  then  fit   empirically  by  an  exponential  series,)   A.e  j 

j=l 

The  reduction  of  the  power  data  to  give  the  appropriate  heat  generation 

rate  in  the  fuel  region  was  accomplished  in  the  following  manner.   The 

t  t'LRT-I  core  contained  28  assemblies,  51  plates  per  assembly  and  an  active 

3 
volume  of  7.523  cm  per  plate.   Therefore, 

H     cm    (P(t)  x  10    .  .    1  assembly..  ,  1  plate  . 

plate       28  assemblies    51  plates    _  roo   3 

r  7 .523  cm 

(16) 

3 
=  97.78  P(t)  watts/cm 

was  the  average  heat  generation  rate  in  an  average  fuel  plate,  where 

P(t)  is  the  total  power  in  megawatts.   Converting  this  to  the  required 

3         —  '3 

dimensions  of  cal/sec-cm  yielded  H  ,    (t)  =  23.37  P(t)  cal/sec  cm  . 

J  plate 

The  heat  generation  rate  in  the  center  of  the  fuel  plate  is  found  in  terms 
of  the  average  heat  generation  rate  since  the  heat  generation  rate  is 
proportional  to  the  neutron  flux  distribution. 
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L  L 

(   C  0  (x)  Adx     f   H  cosh  Kx  dx  H  sinh.  kL 

H     (t)  = =  =  (i7) 

Plate        jL  Adx  L  K  L 


where  L  is  half-thickness  of  plate,  0.0254  cm, 

and  H   is  heat  generation  rate  at  center  of  the 
average  fuel  plate. 


%   was  determined  from  the  neutron  transport  theory  relationship  for 

heavy  absorbers  (44), 

^-  =  tanh-|-  (18) 

tot         s 


to  be  0.7973  cm"   . 


Therefore, 


KL/sinh  KL=  g; 02025  =  1"°  ( 19) 


H  .  (t)  =  H  '    (t)  =  23.37  P(t)  cal/sec  cm3.      (20) 
of        plate 


The  correction  from  the  average  fuel  element  to  the  one  of  interest, 
a  central  fuel  element,  required  a  maximum  to  average  correction.   There- 
fore, 

H  _(t)  =  H  _(t)  (0   /0)  =  23.37  P(t)  (1.9)  =  44.38  P(t)  cal/cm3sec.   (21) 
of       of      max 

The  final  correction  was  to  assume  that  approximately  5%  of  the  power  was 
generated  in  the  moderator.   Therefore, 

H  _(t)  =  0.95  H'  (t)  =  42.16  P(t)  cal/cm3sec.       (22) 
of  of 

The  relation  used  to  calculate  the  heat  generation  rate  in  the  moderator 

(23) 

qm(t)  =  0.05  q£(t)  /=  0.05  q£<B^=  0.05  Hof<t)  (££ff^5>  =  0.017695  Ho£(t) 

m  m 
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where  qf  is  average  heating  rate  in  the  fuel, 

V  is  volume  of  the  fuel, 

V  is  volume  of  the  moderator, 
m 

Lf  is  half-thickness  of  the  fuel, 

L  is  half-thickness  of  the  moderator, 
m 

and    q   is  the  heat  generation  rate  in  the  moderator. 

^m  ° 

The  above  equation  assumes  that  there  is  a  flat  spatial  distribution  of 

the  heating  rate,  that  approximately  5%  of  the  total  heat  generation 

takes  place  in  the  moderator  and  that  the  average  heat  generation  rate 

in  the  fuel  is  well  approximated  by  the  "heat  generation  rate  in  the 

center  of  the  fuel. 

As  previously  mentioned  that  data  from  the  temperature  traces  were 

fit  with  a  finite  number  of  terms  of  a  Fourier  series  of  the  form 

P 


V  k     2nnt 

)   b  cos  

U       n     a 


where 


(24) 
n=0 


bo=f  j*y  (t)  dt, 


b  =  1  a  ja  y(t)  cos  -^  dt, 
n        o  J  a 

y(t)  =  experimental  temperature  trace  data, 

and     a  =  interval  of  periodicity. 

The  data  reduced  from  the  power  traces,  actually  H  f(t),  were  fitted 

with  a  finite  number  of  terms  of  an  exponential  function  of  the  form, 

s 

£  A.  e  V  ,  (25) 

j=l 

where  A,  and  K.  were  parameters  which  were  determined  by  trial  and  error 

to  give  the  best  fit.   The  best  fit  parameters  were  determined  by  means 
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of  an  IBM-6.50  computer  program  described  in  Appendix  B.   This  program 
resulted  from  a  very  minor  modification  of  one  written  by  L.  R.  Foulke  (17) 
The  data  for  H  p(t)  and  the  approximating  equations  are  shown  in  Figures 
8  through  11. 

The  thermal,  nuclear  and  geometric  constants  used  in  determining 
the  temperature  distributions  are  given  in  Table  1. 

Table  1.   Constants  Used  to  Evaluate  the  Temperature  Distributions 


Constant 


Fuel 


Moderator 


a,  Thermal  Diffusivity,  cm  /sec 

K,  Inverse  Diffusion  Length,  cm 

L,  Half-thickness  of  Region,  cm 

k,  Thermal  Conductivity,  cal/cm  sec  C 


0.82  0.001512 

0.7973  0.0 

0.0254  0.071755 

0.5002  0.001488 


2.2.2  Surface  Heat  Flow.   The  heat  flow  rate  out  of  the  fuel  and 
into  the  moderator  as  a  function  of  time  was  evaluated  by  forming  the 
partial  derivative  with  respect  to  position  evaluating  it  at  the  outside 
of  the  respective  region  and  multiplying  by  the  respective  thermal  con- 
ductivity. The  heat  flow  out  of  the  fuel  and  into  the  moderator  are, 
respectively, 


99f(x,t) 

(q/A).(t)  =  -k-  — -jj 

n    f        f   9x 


(26) 


x=L 


and 


(q/A)  (t)  =  -k    m 
m        m  - 


9x 


(27) 


x=L, 
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Evaluating   the  above  equations  yields   in  the   fuel 


(q/A)f(t)   -    -   kf|     ),<y  JL)Bl    (Dt   cos  B±t  +   K1    sin  6^0 


2   2 

n  tt  Q 

(§f)  AL2 

2L     e 


E 


n=l,3,5,-(L     /n^) 


2   2 

P  rn  tt   g,x 


q       a   A.   cosh  kL 
oo  i 


4  4  2 
16L4  i  f     4L2  J        4L2 


2   2  2    2 

njl_a  ^    M&JL&  +  a    2 


+ 


E^ 


x.t 


A4   e     j 


Y     a   cosh(KL)sinh(y  —~ 


L) 


j=l     k.   (a*     -  \.) 
£  j 


:oBh(/-J- 


K    sinh  k 


L) 


«} 


and   in  the  moderator 


(q/A)    (t)  =   -k     /Y   ,    ri^B.    (D,    cos  8,t  +  E.    sin  B.t) 


2   2 

n  tt  g 


CO 


V      (— ) 


4L 
e 


tip  1,31, 5,'"  L  /nna 


E 


a  F  A, 


± 


2    2 

n  tt  g 


2   2 


+ 


j=l      k    (^^r  +   M    C12-^) 


m     4L2 


4L 


2   2 

Bi        4L2 

4  4   2 

n  tt  a        Rz      , 

4  pi 

16L 


-  i  aF  ai   /  vSjw5ll)i  ) 

J=i    *  ^        v      II*  A7       /  J' 


»  j 


V  cosh    /     i 

Ha       } 


where  D,  — 


cosKy  2^  L)sinh(y  ~  L)-  cos(y  ~  L)sln(Y  2a    L) 


cos(y—  Dcosh  (y  — 


.n2(^L)8inhV^L) 
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and  E,   =   COSh^~£  L)slnW  ^  L>  +   CCW^  L^5tnV^  L> 

~^~         "c         T7H7 r"7!T 


cos'(»/- —  L)  cosh  (y  —  L)+-  sin  (y  —  Dsinh  (  y  —  L ) 


2.3  Reactivity  Effects  Due  to  Temperature  Coefficient 
and  Fuel  Expansion 


It  is  pointed  out  by  Deverall  and  Griff ing  (9)  that  the  temperature 
rise  in  the  moderator  in  the  central  unit  cell  cannot  be  used  directly 
to  determine  reactivity  changes.   "Since  the  temperature  coefficient  of 
reactivity,  a,  was  determined  under  conditions  of  a  uniform  temperature 
throughout  the  core  -  a  condition  that  does  not  exist  in  a  transient  - 
it  Is  necessary  to  define  a  properly  weighted  average  temperature.   This 
average  temperature  would  then  produce  the  same  change  in  reactivity  as 
if  an  actual  uniform  temperature  change  of  this  amount  had  been  made. 
This  average  is  defined  by 

2g  =  i  i  (x)  at  d)  dt  (30) 

S    I  (x)  dx 

where 

-+•  ->■ 

AT(x)  is  the  change  of  temperature  at  position  x  , 

■>  -+ 

I(x)  is  the  statistical  importance  at  position  x, 

and  the  integration  is  carried  out  over  the  whole  volume  of  the  reactor." 

The  authors  also  show  that  for  SPERT  I-A  (17/28)  core,  the  system 

under  consideration, 

AT 
TJf—  =  0.65  (31) 

max 

Therefore, 

(32) 
Ak  =  0.65  a(T)  AT 

max 
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A  value  for  a  (t)  of  0.9  x  10"  (Ak/C°)was  used  yielding 


Ak  =(-5.85  x  10"5/  C°)AT    .  (33) 

v  max 


A  similar  problem  was  faced  by  Forbes  (13).   In  determining  the 
reactivity  effect  due  to  fuel  plate  expansion  he  stated, 

"In  order  to  obtain  the  reactivity  change,  the  temperature 
distribution  and' void  importance  function  in  the  core  must  be  combined 
to  obtain  the  dynamic  reactivity  coefficient  as  opposed  to  the  static 
coefficient  which  applies  only  to  uniform  void  distributions.   Applying 
the  observed  distribution  functions  for  temperature  and  void  worth,  it 
is  found  that  the  effective  average  temperature  rise  under  dynamic 
conditions  can  be  obtained  from  the  temperature  rise  at  the  center  ;of 
the  core  by  the  relation 

A0  =  0.7  AG    .  (34) 

max 

The  reactivity  change  due  to  plate  expansion,  Ak.,will  be 

Ak,  =(p)  3  a  v  (0.7  A3   ),  (35) 

1   ov  max 

3k 
where (r — )  is  the  average  void  coefficient  for  the  core, 

a  is  the  linear  expansion  coefficient  of  aluminum, 

v  is  the  volume  of  the  aluminum  which  is  heated  (i.e.,  the 
volume  of  fuel  plates  proper), 

AG"  is  the  average  temperature  rise  of  the  aluminum 

and  AG    is  the  temperature  rise  at  center  of  the  core. 

max 

For  the  SPERT  I-A  (17/28)  core  the  appropriate  constants  are  the  following: 
a  =  2.5  x  10"5/°C 
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4   3 
v  =  2.8  x  10  cm 

(|^)=  -3.5  x  10"6  Ak/cm3 
ov 

Therefore,  the  expression  for  the  reactivity  change  becomes 

Ak,  =  (-5  x  10"6  ^  )  (AG   ).  "  •  (36) 

1  o       max 

In  order  to  avoid  erroneously  taking  into  account  the  void  formation 
due  to  fuel  element  expansion  in  both  the  temperature  coefficient  and 
in  the  fuel  element  expansion  calculation,  the  fuel  element  expansion 
was  calculated  only  for  the  temperature  rise  in  the  fuel  over  the 
temperature  rise  in  the  moderator.  Therefore,  the  reactivity  effects 
due  to  the  temperature  coefficient,  Akp,  and  due  to  the  fuel  plate 
expansion,  Ak„,  are 

Akj,(t)  =  -5.85  x  10"5  (Ak/°C)  0    (t)  (37) 

and       AkgCt),  =  -5  x  10"6  (Ak/°C)  (6f   ^t)  -  0mod  (t)  )      (38) 


where     0   ,  (t)  is  the  average  temperature  rise  in  the  moderator  at 

time  t 

0_   T(t)  is  the  average  temperature  rise  in  the  fuel  at  time  _t 
as  obtained  from  Tables  2  through  5. 

2.4  Reactivity  Effects  Due  to  Steam  Formation 

The  calculation  of  the  steam  production  was  based  on  the  same 
assumptions  as  those  used  in  the  "Conduction  Boiling  Model  for  Reactor 
Self -Shutdown"  suggested  by  S.  G.  Forbes  (15).   In  Forbes'  work  the  steam 
volume,  V  ,  was  assumed  to  be  proportional  to  a  fraction,  f   ,  of  the  energy, 
E  ,  transferred  to  the  moderator  after  the  time  boiling  first  occurred 
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in  the  core.  The  steam  volume  is  given  by 

X,     -  £    E 

Vs  -  _*L_JL-   f  (39) 

8 

where  h  is  the  energy  required  to  form  a  unit  volume  of  steam  from 
s 

3 
boiling  water  at  standard  pressure  (1.35  watt-sec/cm  steam).  The  term 

f   was  regarded  as  a  combination  of  factors  involving  the  fraction  of 

the  energy  actual'ly  forming  steam  during  nucleate  boiling  (about  1%) 

and  the  fraction  of  the  core  heat  transfer  area,  A,  which  is  involved 

in  boiling  heat  transfer  (about  10%)  (15).  The  reactivity  effect  of 

the  steam  is 

Ak  =  V  C  =  fas  gs  r  ,._. 

s    s  v   — r C   ,  (40) 

h     v 
s 

where  C   is  the  void  coefficient  in  the  center  of  the  core,  (C  -'1.2   x 
v  v 

-6      3 
10  Ak/cm  for  Spert  I-A) .   In  this  investigation  the  factor  f   was 

3.  S 

divided  into  its  two  components,  the  fraction  of  the  core  involved  in 
boiling,  f  ,  and  the  fraction  of  the  energy  actually  forming  steam  during 
nucleate  boiling,  f  .   This  was  done  since  it  was  possible  to  approximate 

S 

f  directly  from  the  fuel  surface  temperature  and  the  assumption  that 

the  gross  temperature  distribution  over  the  core  was  proportional  to  the 

bare  core  power  distribution.   It  is  expected  that  the  final  factor,  f  , 

will  be  independent  of  the  pulse  parameters  for  a  particular  system  and 

that  it  will  prove  essentially  independent  of  the  reactor  parameters  in 

any  heterogeneous  water  moderated  system.   In  effect,  the  final  parameter, 

f  ,  was  left  to  be  calibrated  by  any  particular  pulse.  The  test  of  this 
s 

model  was  of  course  a  constant  f  .   For  the  two  boiling  runs  considered, 

-3  -3 

the  values  of  f  calculated  were  3.7  x  10  '  and  3.3  x  10  "  which  differ 
s 
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by  less  than  270.  The  27o  difference  is  less  than  would  be  expected  for 
the  accuracy  of  the  Input  data.  The  reactivity  effect  of  the  steam  is 
then 


(AD 


Ak  (t)  =(~7,28  X  10  ,  Ak/MW-sec)  (f  f   E  (t)  ) 
1.35  x  ID"6  a  S   S 


=  -1.87  x  10   f  E  (t) 
a  s 


f  Is  fractidn  of  core  heat  transfer  area  involved  in  boiling, 

f  is  fraction  of  energy  actually  forming  steam, 
s 

and  E  (t)  is  total  energy  into  core  after  initial  boiling  in  MW-sec. 

s  / 

The  fraction  f  was  calculated  assuming  the  gross  core  temperature 

cl 

distribution  had  reached  a  dynamic  equilibrium  with  the  power  and  the 
power  distribution  could  be  approximated  by  that  of  an  equivalent  bare 
core  with  an  effective  height,  2Z  ,  and  radius,  R  .   The  gross  temperature 
distribution  in  the  core  is  then 


Q  ,   „>    n   T  ,2.A0A8  rN       ttZ 
0  (r,Z)  =  0O.JO(— )   cos  — 


(A2) 


where  0o-is  the  surface  temperature  at  the  time  of  interest  on  the  axial 

center  line  of  a  central  fuel  element.   The  maximum  value  of  Z,  Z 

max 

for  which  boiling  will  occur  on  any  plate  can  then  be  obtained  knowing 

the  boiling   temperature,  0,  ,  and  the  axial  centerline  surface  temperature, 

e  J(^A8r). 
o  o   R 


That  is 


ttZ 


cos 


max 


2  Z 


0   J  (240A8r} 


o  o   R 


(A3) 


33 


so  that 

Z  o  i  9u 

max  _  £       -1  , b v 

Z       tt    C°8    a   T  ,2.4048r, 

e  0   J  ( — - ) 

o  o   R 
e 

The  maximum  value  of  r  for  which  boiling  will  occur,  r   ,  on 

—  °  max 

the  core  axial  centerline  can  be  obtained  from 

a  _  a    r   2.4048r 

0,  =  U   J   , max  x 

boo  ( r ) 

e 

so  that  (44) 

max  _    1     J  r_b* 

R     2.4048    °{Q 
e  o 

where  y  =  J    (x)  is  the  inverse  of  x  =  J  (y) . 
J  o  o  J 

The  volume  fraction  of  the  core  having  a  fuel  surface  temperature 
above  the  boiling  temperature  is 

1  r 

f     =        D2   m7   x      5   maX      2  Z  2rrrdr 

v  ttR^   ( 2Z   )      Jo  max 

e  e 

=  _g_        'ma*         z^ 

R   2        °  Z—  rdr 

e  e 

Making  the  change  in  variable  n,  =  r'R  ,  the  integral  becomes 

n      Z   (n) 

f  =2(  maX    T      ndn  .  (46) 

v     o      z 
e 

Substituting  the  value  of  Z    /  Z  from  equation  (43)  yields 
b  max  '   e 

n  ,  b 

max  -1 


f     =   4/tt   5   max     cos"1     Q      T  ,2.4048   ,        ndn    .  (47) 

v  Jo  0     J  (  n)        '    ' 

oo    o 

This  integration  was  then  carried  out  numerically  using  the  value 

of  n  ,   -  r   /  R  determined  from  equation  (44).   Since  there  is  a 
max    max   e 
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constant  heat  transfer  area  per  unit  volume  in  the  core  then  f  -  f  . 

_v   __a 

For  the  two  boiling  runs  t=  15.8  msec  and  23  msec,  f  was  equal  to 
0.116  and  0.084,  respectively. 

The  total  energy  transferred  to  the  moderator  after  the  time  of 
initial  boiling  was  obtained  by  considering  the  moderator  volume 
associated  with  unit  surface  area  in  the  central  fuel  element.   The  heat 
content  of  the  mpderator  at  the  time  boiling  temperatures  were  reached 
at  the  surface  and  at  the  time  of  peak  power  were  calculated  based  on 
the  conduction  model.   While  this  model  gave  a  somewhat  erroneous 
temperature  distribution  above  boiling  temperatures  it  accurately  repre- 
sented the  heat  flow  into  the  moderator.   The  difference  between  the 
moderator  heat  content  at  the  time  of  interest  and  at  the  time  boiling 
temperatures  were  first  reached  was  the  energy  available  for  boiling  per 
unit  fuel  surface  area.   Plots  of  the  central  fuel  surface  temperature 
and  the  average  moderator  temperature  in  a  central  element  used  to 
calculate  the  moderator  heat  contents  are  shown  in  Figures  2  and  3.   The 
total  energy  available  for  steam  formation  is  obtained  by  multiplying 
by  the  total  heat  transfer  surface  area  of  the  core.   This  somewhat  over- 
estimates the  total  energy  but  is  probably  the  best  estimate  of  the 
energy  of  interest  since  the  boiling  region  is  confined  to  a  rather  small 
central  portion  of  the  core.   Therefore 

E  (t)  =  (9  (t)  -*9  (tu))  C  M 
s        m       m  b    p 

where  C   is  the  heat  capacity  of  the  moderator 
P 
M  is  the  mass  of  the  moderator  in  the  system 

9  (t)  is  the  average  moderator 
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Figure  2.   Graph  used  to  determine  the  moderator  energy 
for  boiling  calculations  during  a  transient  with 
an  initial  period  of  95.8  msec. 
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Figure  3.    Graph  used   to  determine  the  moderator  energy    content 
for  boiling    calculations  during  a  transient  with  an  initial 
period    of    23  msec. 
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3.0  RESULTS  AND  DISCUSSION 
3.1  Temperature  Distribution  and  Surface  Heat  Flow 

The  temperature  distributions  as  obtained  from  equations  15  and 
16  are  shown  in  Figures  A  through  7.   These  plots  were  obtained  using 
experimental  data  from  four  transient  tests  on  the  SPERT  I-A  (17/28) 
reactor  having  initial  periods,  t,  of  15.8,  23,  120,  and  150  msec, 
respectively.   The  transient  burst  for  t  equal  15.8  and  23  msec 
show  regions  In  the  moderator  which  have  temperatures  above  the  boiling 
point  at  a  pressure  of  one  atmosphere.  .It  is  not  believed  that  this 
superheating  takes  place.   These  temperature  distributions  are  shown 
since  such  a  small  portion  of  the  moderator  is  above  the  saturation 
point  that  it  is  not  likely  that  it  will  materially  affect  the  temperature 
in  the  remainder  of  the  moderator  or  the  average  moderator  temperature. 
One  possibility  which  must  be  considered  in  calculating  the  reactivity 
effects  If  that  pressure  transients  are  developed  which  raise  the  boiling 
point  above  the  temperatures  observed  in  the  core.   This  appears  not  to 
be  the  case  for  two  reasons.   First,  experimental  measurements  of  the 
pressures  do  not  indicate  sufficient  rises  in  pressure  and  second,  the 
total  reactivity  compensations  at  peak  power  indicate  that  boiling  must 
have  taken  place. 

The  experimental  surface  temperature  traces  and  the  approximate 

analytical  fits  for  the  four  transient  tests  are  shown  in  Figures  8  through 

P 

11.   The  parameters  for  the  analytical  fits,  0(L,t)  =  )     B   cos  P  t  , 

1=1 
are  shown  in  Table  2.  The  experimental  power  traces,  actually  H  f(t), 

and  the  approximate  analytical  fits  for  the  same  four  transient  tests 
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Figure  5.   Temperature    distributions,  0(x,t)f  vs  position  in  fuel  and  moderator 
based   on  pure   conduction  during  a   transient  with  an  initial  period 
of    23  msec. 
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Table    2.      Numerical   Values   of   Parameters   for  Emnirical   Fits 
of   A(L,t)    Used    in   Equations    (15)    and    (16) 


q;,msec 


150 


14.593 


o.on 


-20.. 958 


3.1416 


+9.474 


6.2832 


3.763 


9.4279 


+1.313 


12.566 


-0.871 


15.708 


f0.582 


18.850 


-0.447 


21.991 


120 


+12.033 


0.00 


15.831 


^4.760 


+5.875 


9.520 


-2.067 


.4.280 


+6.817 


19.040 


-0.829 


23.800 


f0.399 


'28.560 


■0,100 


33.320 


23 


15.  8W 


+35.752 


o.oo 


42.074 


0.00 


54,937 


17.405 


•59.957 


34.906 


+26.246 


34.810 


+25.821 


69.812 


■7.248 


2.118 


+4.476 


2.215 


-9. 279 


104.718 


+1.094 


139.624 


37.024 


0.838 


174.530 


-3.802 


104.429 


+  2.413 


121.834 
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are  shown  in  Figures  12  thro?  gh  15.  The  parameters  for  the  analytical 

fits,  H  c(t)  =  )   A,   J   are  shown  in  Table  3. 
of      L       j 

One  physical  check  on  the  solutions  not  required  by  the  mathematical 
formulation  of  the  problem  is  that  the  heat  flow  out  of  the  fuel  must 
equal  the  heat  flow  into  the  moderator.  The  heat  flow  data  are  shown  in 
Figures  16  through  19.   It  is  obvious  from  inspection  of  these  data  that 
the  equivalent  heat  flow  condition  is  not  well  satisfied.  There  are 
several  possible  explanations  for  this  discrepancy.  First,  the  heat 
flow  data  is  somewhat  more  sensitive  to -the  analytical  fits  of  the  sur- 
face temperature  and  the  heat  generation  rates  than  the  average  temperatures. 
Second,  the  cladding  between  the  meat  and  the  moderator  was  neglected  in 
calculating  the  temperature  distributions.   Again,  the  average  tempera- 
tures are  far  less  sensitive  to  this  approximation  than  the  heat  flow 
calculations.   Finally,  it  can  be  seen  by  investigating  the  heat  flow 
equations  that  the  discrepancies  could  be  decreased  by  introducing  a 
positive  phase  angle  to  the  surface  temperature  fits.   This  could  be 
attributed  to  a  delay  time  in  the  surface  temperature  measurements. 

The  average  temperatures  in  the  fuel  and  moderator  as  a  function 
of  time  were  calculated  by  means  of  a  numerical  integration  of  the  cal- 
culated temperature  distributions.   These  data  are  given  in  Tables  A 
through  7 . 

3.2  Reactivity  Effects 

The  reactivity  compensations 

Ak  (t)  =  Ak(o)  -  Ak(t) 
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Table   3.      Numerical  Values  of  Parameters   for  Emoirical  Fits 

of  H    -(t)    Used   in   Equations    (15)    and   (16) 
of 


"t;,msec 


150 


+1.912 


7.426 


2.097x10' 


21.090 


+6.062x10"** 


;27 .  260 


120 


9.218 


'9 .  204 


■6.486x10' 


20.233 


+3.564x10 


30.093 


-1.856x10 


42.208 


+3.019x10 


T9 


85.147 


23 


2.414x10 


1 


39.761 


-4.364x10 


124.03 


9.563x10 


10 


187.61 


15.8 


■2.993 


104.11 
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t ,  sec 


Table  4.   Average  Temoerature  Rise  In  the 

Fuel  and  Moderator  for  %  -    15.8  msec. 


t , sec. 

3,      At)   °C 
fuel 

S    /t),oC 

mod 

("§.    n-S     ,),°c 

fuel     mod 

.085 

138.53 

16.732 

121.80 

.080 

127.34 

13.99 

113.35 

0.075 

110.34 

10.81 

99.53 

.070 

89.99 

8.664 

81.33 

.060 

50.12 

4.87 

45.25 

.050 

24.58 

2.446 

22.13 

.040 

13.91 

1.333 

12.58 

.030 

8.51 

0.692 

7.82 

.020 

3.52 

.271 

3.25 

.010 

1.039 

.  .0767 

.962 

Table  5.   Average  Temoerature  Rise  in  the  Fuel 
and  Moderator  for  i  -    23  msec. 


,    At),°c 

fuel 

e    ,(t),  °c 

mod 

(9.    ,-e    ,)  °c 

fuel      mod    , 

117.76 

16.832 

100.93 

108.71 

13.371 

95.34 

86.92 

9.764 

77.16 

59.48 

6.459 

53.02 

36.58 

4.218 

32.35 

23.17 

2.759 

20.41 

16.28 

1.869 

15.21 

11.25 

1.231 

10.02 

3.13 

.260 

2.87 

.16 
.15 
.14 
.13 
.12 
.11 
.10 
.09 
.05 
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Table  6.   Average  Temperature  Rise  in  the  Fuel 
and  Moderator  for  <r  =  120  msec. 


0£     ,(t),°c 

9       ,(t),°C 

sec . 

fuel 

mod 

.59 

32.66 

8.572 

.57 

29.95 

7.738 

.55 

27.33 

6.935 

.50 

21.15 

5.212 

.45 

16.40 

3.834 

.40 

11.98 

-      2.660 

.35 

7.84 

1.744 

.30 

4.89 

1.143 

.20 

2.59 

.534 

.10 

1.17 

.164 

(9.    ,-e    ,),°c 

fuel   mod 


24.09 

22.21 

20.39 

15.94 

12.57 

9.32 

6.10 

3.75 

2.06 

1.01 


Table  7.   Average  Temperature  Rise  in  the  Fuel 
and  Moderator  for  x  =  150  msec. 


sec . 

e.    -(t),°c 

fuel 

6       ,(t),°C 

mod 

.90 

44.46 

12.994 

.85 

37.49 

10.630 

.80 

30.87 

8.510 

.75 

25.23 

6.678 

.70 

20.31 

5.111 

.65 

15.70 

4.644 

.60 

11.48 

2.750 

.55 

8.11 

1.989 

.50 

4.86 

1.472 

.40 

3.58 

.834 

.30 

1.91 

.431 

.10 

.602 

.0652 

(0: 

fui 

=  1 

)       ,),°C 
mod 

31, 

,47 

26, 

,86 

22, 

,36 

18, 

,55 

15, 

,20 

11, 

,06 

8, 

,73 

6, 

,12 

4, 

,39 

2, 

,75 

1, 

,48 

,537 
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where  Ak(o) ■ is  the  initial  reactivity  insertion  to  start  the  transient 

and  Ak(t)  is  the  excess  reactivity  of  the  system  at  any  time,  J:,  are 

shown  as  a  function  of  time  in  Tables  8  through  11.   Tables  8  through 

11  also  show  the  components  of  the  compensated  reactivity  due  to  the 

temperature  coefficient  Ak^tt),  fuel  plate  expansion  Ak  (t)  and  steam 

formation,  Ak  (t).   The  excess  reactivity  Ak( t)  for  two  of  the  transients, 
s 

T  equals  120  and  150  msec  respectively,  is  compared  with  equivalent  data 
obtained  from  a  kinetic  analysis  of  the  power  burst  shapes  by  Miller 
(34)  in  Figures  20  and  21.   The  reactivity  compensation  at  peak  power, 
is  shown  along  with  comparable  data  from* the  kinetic  analysis  in  Figure  22. 
Figure  22  also  includes  the  components  of  the  reactivity  compensations 
for  a  model  suggested  by  S.  G.  Forbes  (15)  as  well  as  for  the  model 
suggested  in  this  report. 

3.3  Conclusions 

The  forms  of  all  of  the  solutions  shown  in  equations  (9)  through 
(16)  are  such  that  three  terms  are  developed.   The  first  term  represents 
the  steady  state  solution  resulting  from  the  surface  temperature  boundary 
condition.   The  second  term  includes  the  transient  portion  of  both  the 
surface  temperature  boundary  condition  and  the  forcing  function,  the 
heat  generation  rate.  The  final  term  represents  the  steady  state  or 
equilibrium  solution  resulting  from  the  forcing  function.   It  has  been 
common  in  several  previous  works  (27,36)  to  assume  that  the  temperature 
is  separable  in  space  and  time.   It  can  be  seen  from  the  derived  solu- 
tions that  this  will  be  a  good  approximation  of  the  temperature  distribution 
only  when  the  second  term,  the  transient  solution,  has  died  out.   The 
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Figure  20.  Comparison  of  calculated   percent  excess 
reactivities,  Ak,  and  those  obtained  by  kinetic  analy- 
sis  vs  arbitrary   time  during  a  transient  with  an 
initial  period  of  120  msec. 


0.72 

070 

0.68 

^0  66 

• 

^n 

T™ 

< 

r~~ — ' 

^^ 

^"V^i 

> 

> 

pi 

^0.64 
062 

\.    < 

e 

—            < 

0.60 

—  KINETIC    ANALYSIS     IDC 
»     CALCULATED    VALUES 

1    1    1    1    1    1 

(-16317 

( 

) 

ti 

0 

.1 

0 

2 

0 
TIM 

3 
E,SE 

0 

•a 

.4 

0 

5 

0 

6 

Figure  21.  Comparison  of  calculated  percent  excess 
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Table  8.   Compensated  Reactivities  for 
T7  =  15.8  msec  Run 


t  ,sec 

AkpCt) 

Akg(t) 

Akg(t) 

Akc(t) 

0.084 

0.979  x 

io"3 

0.609  x   10~3 

1.86  x 

io"3 

3.45   x 

io"3 

0.080 

0.818  x 

10"S 

0.546  x   IO"3 

1.21   x 

io"3 

2.60  x 

io"3 

0.075 

0.632  x 

io"3 

0.498  x   IO-3 

0.46  x 

io"3 

1.59  x 

io"3 

0.070 

0.684  x 

io"3 

0.407  x   IO"3 

0.13  x 

10   3 

1.22  x 

io"3 

0.060 

0.285  x 

io"3 

0.226  x    10~3 

0.0 

0.51   x 

io"3 

0.050 

0.143   x 

io"3 

0.111  x   IO-3 

0.0 

0.25  x 

io"3 

0.040 

0.078  x 

io"3 

0.063  x   10~3 

0.0 

0.14  x 

io"3 

0.030 

0.050  x 

io"3 

0.039  xio"3 

0.0 

0.08  x 

io-3 

0.020 

0.016  x 

io"3 

0.016  x   10~3 

0.0 

0.03  x 

io"3 

0.010 

0.004  x 

10~3 

0.005  x   io"*3 

0.0 

0.01   x 

io-3 

Table  9.   Compensated  Reactivites  for 
%  =  23  msec  Run 


t ,  sec 

Ayt) 

AkE(t) 

Aks(t) 

Akc(t) 

0.16 

0.985  x 

io"3 

0.505   x 

io"3 

1.21   x   10"3 

2.70  x   10~3 

0.15 

0.782  x 

io"3 

0.468  x 

io-3 

0.54  x    10~3 

1.79  x   IO"3 

0.14 

0.571   x 

io"3 

0.386  x 

io-3 

0.003  x   IO"3 

0.96   x   10  3 

0.13 

0.378  x 

io"3 

0.265  x 

10  3 

0.0 

0.64  x   10~3 

0.12 

0.247  x 

10   3 

0.162  x 

10   3 

0.0 

0.41   x    10~3 

0.11 

0.161   x 

io"3 

0.102  x 

10~3 

0.0 

0.26  x   IO-3 

0.10 

0.109  x 

io"3 

0.076  x 

io"3 

0.0 

0.18  x   IO"3 

0.09    " 

0.072  x 

io-3 

0.050  x 

io"3 

0.0 

0.12  x   IO-3 

0.05 

0.015  x 

io"3 

0.014  x 

io"3 

0.0 

0.03  x   IO"3 

63 


Table  10.   Compensated  Reactivities  for 
1  -    120  msec  Run 


t  ,sec 

Aiyt) 

Akj(t) 

Akg(t) 

Akc. 

:t) 

1 

0.59 

0.501  x   10"3 

0.120  x 

io"3 

0.0 

0.621 

X 

io'3 

0.57 

0.453   x   10"3 

0.111   x 

10"3 

0.0 

0.564 

X 

10~3 

0.55 

0.406  x   10~3 

0.100   x 

io'3 

0.0 

0.506 

X 

io"3 

0.50 

0.305  x    10~3 

0.080  x 

io"3 

0.0 

0.385 

X 

io"3 

0.45 

0.224  x   10"3 

0.063  x 

IO"3 

0.0 

0.287 

X 

10~3 

0.40 

0.156  x    10"3 

0.047   x 

IO"3 

0.0 

0.203 

X 

IO"3 

0.35 

0.102  x   10-3 

0.030  x 

io"3 

0.0 

0.132 

X 

io"3 

0.30 

0.067  x    10"3 

0.019  x 

io"3 

0.0 

0.086 

X 

io"3 

0.20 

0.031   x   10"3 

0.010  x 

io"3 

0.0 

0.041 

X 

io"3 

0.10 

0.010  x    10"3 

0.005  x 

IO**3 

0.0 

0.015 

X 

io"3 

Table  11.   Compensated  Reactivities  for 
%  -    150  msec  Run 


t  ,sec 

AkT(t) 

Ak  (t) 

i 

Akg(t) 

Akc(t) 

0.90 

0.760 

X 

io"3 

0.157  x 

io"3 

0.0 

0.915   x    IO"3 

0.85 

0.622 

X 

io"3 

0.134  x 

io"3 

0.0 

0.756   x    IO-3 

0.80 

0.498 

X 

io"3 

0.112  x 

io"3 

0.0 

0.610  x   10~3 

0.75 

0.391 

X 

io"3 

0.093  x 

io"3 

0.0 

0.484  x   IO"3 

0.70 

0.299 

X 

io"3 

0.076  x 

10~3 

0.0 

0.375  x    IO-3 

0.65 

0.272 

X 

io"3 

0.055  x 

io"3 

0.0 

0.327  x   IO-3 

0.60 

0.161 

X 

io"3 

0.044  x 

io"3 

0.0 

0.205  x    IO"3 

0.55 

0.116 

X 

10~3 

0.031   x 

io"3 

0.0 

0.147  x   IO'3 

0.50 

0.086 

X 

io"3 

0.022  x 

10  3 

0.0 

0.108  x   IO-3 

0.40 

0.049 

X 

io"3 

0.014  x 

io'3 

0.0 

0.063  x    IO"3 

0.30 

0.025 

X 

10~3 

0.007  x 

io"3 

0.0 

0.032  x   IO"3 

0.10 

0.004 

X 

io"3 

0.003  x 

10  3 

0.0 

0.007  x   IO-3 
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numerical  results  obtained  from  evaluating  the  set  of  equations  (9)  and 
(10)  and  the  set  of  equations  (15)  and  (16)  show  that  the  transient  term 
is  negligible  for  all  times  of  interest  in  the  fuel  region  b ut  it  makes 
a  significant  contribution  for  all  times  of  interest  in  the  moderator. 

In  one  sense  it  would  be  more  informative  to  investigate  the  reactor 
burst  behavior  using  the  minimum  of  input  data  (i.e.  the  physical 
dimensions  and  the  initial  reactivity  insertion)  and  test  for  corroboration 
of  all  of  the  experimentally  measured  variables.   However,  it  seemed 
better  in  analyzing  for  the  shutdown  mechanisms  to  use  as  much  of  the 
data  as  possible'  leaving  only  the  compensated  reactivities  as  a  check 
on  the  validity  of  the  model.   The  compensated  reactivity  was  established 
as  a  criterion  because  of  its  extreme  sensitivity  to  the  state  of  the 
system  and  because  of  its  direct  influence  on  the  safety  of  nuclear  reactors, 

The  results  of  this  work  are  two-fold.   First,  a  more  accurate  view 
of  distribution  of  the  energy  during  a  transient  burst  is  presented  and 
second,  "a  model  based  on  the  energy  distribution  was  shown  to  predict 
the  reactivity  effects  as  well  as  any  of  the  existing  models.   The 
advantage  of  this  model  is  that  assuming  the  fractional  energy  associated 
with  void  formation  can  be  determined  as  the  mechanism  of  transient  boiling 
becomes  better  understood,  the  final  empiricism  can  be  removed  from  the 
model. 

3.4  Further  Investigation 

There  are  several  avenues  of  attack  for  further  work  in  determining 
the  inherent  shutdown  mechanisms.   First,  additional  data  on  Spert  I-A 
should  be  tested  on  the  model  proposed  in  this  report  to  make  certain 
that  it  is  as  capable  of  determining  reactivity  effects  as  these  prelim- 
inary runs  indicate.   Second,  application  of  this  model  to  any  new  system 
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will  mean  that  the  surface  temperatures  and  power  traces  would  not  be 
available.   This  problem  can  be  circumvented  by  studying  the  two  region 
conduction  problem  subject  only  to  the  heat  generation  rate  forcing 
function.   The  heat  generation  rate  can  be  calculated  from  the  initial 
reactivity  insertion,  allowing  a  feedback  from  the  induced  negative 
reactivity  to  the  heat  generation  rate  through  the  reactor  kinetics 
equations.   This  , suggests  an  analog  solution  or  possibly  a  digital 
analog  combination.   Third,  application  of  the  heat  transfer  equations 
developed  in  this  report  should  be  used  to  determine  the  mode  of  heat 
transfer  during  transient  operation  by  investigating  a  single  plate  ir 
as  much  detail  as  possible.   The  application  of  this  study  can  probably 
be  done  more  simply  using  electrical  heating.   Fourth,  a  detailed  study 
of  nucleate  boiling  at  low  heat  fluxes  is  necessary  before  complete 
understanding  of  the  mechanisms  of  shutdown  can  be  obtained.   This 
study  should  provide  a  direct  measurement  of  the  fraction  of  the  energy 
used  to  produce  steam,  f  .   Fifth,  work  on  this  model  should  be  extended 
to,  investigate  further  available  evidence  on  other  Spert  reactors  to  see 
if  it  will  account  for  changes  in  other  parameters  such  as  neutron  life- 
time, pressure  and  coolant  flow.   Finally,  experimental  and  analytical 
work  should  be  done  on  the  zirconium  hydride  moderated  Triga  systems 
since  qualitatively  they  show  the  greatest  inherent  safety  that  has 
been  demonstrated  to  date. 
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APPENDIX  A 

Derivations  of  Solutions  for  the 
Temperature  Distribution  in  the 
Fuel  and  Moderator  of  a  Unit  Cell 


The  differential  equations  governing  the  time  dependent  temperature 

distribution  in  the  fuel  and  moderator  are 
s 


^2.  ,   _.   V  q  cosh(rx)a.e  j      ,    „A/   .» 
V  0  (x,t)  +  ).  Jpo i   J   _  j.   aQ(x,t) 

f      .  ,        k        "  a     8t 


(A-l) 


and 


L  V   FA,  eY    ,  a0mu,t) 


2A  /    n   V    F  A.  e  j 


j=l  a    9t 


(A-2) 


The  differential  equations,  (A-l)  and  (A-2),  are  most  easily 
solved  by  use  of  Laplace  transforms.   Considering  the  fuel  region  first 
and  transforming  the  time  variable  in  equation  (1)  yields 

8  ' 
2  a  /    s  ,  V       q   co8h(Kx)  A.    . 
r,  $f(x,s)+^  %o UI[S^(K,S)  -9  U,o)].   (A-3) 

j=i   k  (s  -  \.)    a    £ 

0(x,S j  ie   the  transform  of  9(x,t) ,  and  9(x,o)  is  the  initial  temperature 
distribution.   Assuming  the  initial  temperature  distribution  to  be  flat, 
the  equation  can  be  normalized  by  letting  Q(x,o)  equal  zero.   Thus 
Q(x, t)  is  the  temperature  excess  in  the  fuel  over  the  initial  temperature, 
The  initially  flat  temperature  distribution  in  the  fuel  is  unreal  but  it 
is  a  good  approximation  if  the  transients  are  started  from  low  power 
levels. 

Letting  9(x,o)  =  0  ant  rearranging  the  terms  in  equation  (A-3) 
yields 


V  2  0f  (x,3)    -  |  t>t    (Xi  S)  =   -  £       %oCosh(l<x)A1 

j=l  k    (S    -    \   ) 


(A-4) 
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2 
The  usual  one-dimensional  slab  assumption  is  made  and  V   becomes 

g2 

— -   .   This  assumption  neglects  the  axial  flow  of  heat  in  the  fuel  corn- 
ax2 
pared  to  the  radial.   While  this  might  be  the  limiting  assumption  in  the 

analysis,  it  is  probably  not  seriously  in  error. 

The  homogeneous  solution  to  the  one  dimensional  form  of  equation 

(A-4)  is  well  known  and  is  derived  in  many  standard  texts.   This  solution 

is 

0-,  (x,S)  =  A  cosh  /-  x  +  B  sinh  /-  x  (A-5) 

fh  si  a  n/  a 

The  particular  solution  is  found  by  the  method  of  undetermined  co- 
efficients and  is 

<-,   q   a  A.  cosh  k  x 
<S>-      (x,S)  =  )   -22 3 (A-6) 

tP  .V    k(S-X.)(S-QK  ) 

j  - 1     j 

Therefore, 

r~,  r^,        „   q   a  A.  cosh  fK  x) 

$.   (x,S)  =  A  cosh  /-  x  +  B  sinh  /-  x +)     °°   ,Jrg : .  (A-7) 

fp  v  a  n/  a    L>      k(S-\.)  (S-aiO 


J=l 


wj       2' 


The  boundary  conditions  used  to  determine  the  constants  A  and  B 
are  as  follows   First,  the  temperature  gradient  in  the  center  of  the 
region  (x  =  0)  is  zero  for  all  time.   Second,  the  surface  temperature 

is  matched  with  the  experimental  data.   The  surface  temperature  is 

f1     R  t 
expressed  as  a  sum  of  exponentials,   )   B  e'  i   or  a  Fourier  series, 

)      B  cos  p.t.   Firsr  consider  the  solution  in  which  the  exponential 

i=l 
boundary  condition  is  used. 

The  transformed  boundary  conditions  are 
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X{^o    }<{o 


e,(x,t         .^u  or    dif(x,s^  ^  (A_8) 


9x  '  dx 

x=0  x=0 


and 


P 


^{ef(L,t)}     =  £{[Bl   JV}     or     »f<lt»=£ 


B 


i  .       (A-9) 


i=l  i^l(S-Pi) 

Applying  eq.uation  (A-8)  to  equation  (A-7)  shows  that  B  must  equal 
zero.   A  is  determined  by  evaluating  equation  (A-9).   The  solution  in 
the  transform  domain  is  shovn  in  equation  (A-10). 

-    ,     „N        f1  B      cosh     /-  x  i,         q       aA.  r 

Vx's)  =i         -i *-* -  £      -^_ J _   {  cosh  (KX)  . 


•-i  ,c   □    ^  u      /^\  k(S-\.) 

i=l  (S-F$.)    cosh     /-  L      ._,  1 

■  i  \'  a  j-1 


(S-aK  ) 


S 
cosh    (kL)      coshv/_Q_x      1  (A-10) 


cosh  y  a 


I 


Transforming  equation  (A-10)  back  into  the  time  domain  is  greatly 
simolified  since  no  poles  of  order  greater  than  one  occur  in  the 
inversion  integral.   It  is  interesting  to  note  that  it  is  the  occurence 
of  higher  order  poles  which  conmlicates  the  solution  to  the  multiple 
region  nroblem.   Equation  ( A- 11) , shovn  below,  (33)  can  be  used  to  invert 
transformed  functions  of  the  form  f(s)  =  j(s)  /  j^(s)  if  the  degree  of 
J^(s)  is  at  least  one  greater  than  j(s)  and  only  poles  of  order  one 

occur, 

m 

ifn    ^ 

(A-ll) 


/    -1  jfU)}     =    Y         ""n'  P   t 
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p   denotes,  the  n  simple  poles  of  f  (  S)and  Jl(^n)    denotes  the  valup  of 

dl(S)  _ 

— ttt-  evaluated  at  »=  p  . 
dS  _n 

The  inversion  of  equation  (A- 10)  is  shown  in  detail  below.   First 
consider  the  third  term  of  equation  (A-10) . 

f-\{  y     qroa  cosh  KL  Aj    cosh  ^x  s     q^   fl  cosh  (  kL)  A 

*    \-i  2  /F~  =  L  E ^'    ] 

J-1  k(s-K  )(S-aK  )  cosh  /-  L     j=l 

2 
The  obvious  poles  of  I  are  at  S  =  X.  and  S  =  gx      .   Additional 

roots  of  the  denominator  exist  for  cosh  /—  L  =  0.   The  roots  of 

,v  a 

cosh  /—  L  are  obtained  by  making  the  transformation  S  =  -  yd* .      This  is 

v  a 

done  since  cosh  /—  L  cannot  equal  zero  for  real  values  of  the  argument 
and  a  and  L  are  always  positive,  real  constants. 


(A-12) 


/T7  ZZT        /GT 

cosh  /—  L  =  cosh  i  /—  L  =   cos  /   L  =  0 
v  a         v  xa        v  a 


(A-13) 


therefore 


—  R  -  —  ♦  n  odd 
a  o    2 

2  2 

j  _  "  n  q 

*n   "   /T2  (A- 14) 

4L 


2  2 
4L2 


Evaluating^    -I I  >   by  means  of  equation  (A- 11)  at  the  poles 

^     2  2 

p  =  gv      ,  p  =  K   and  P  =  -  ^  yields  equation  (A-15). 

4L 


76 


y-l  r  \       _  cosH-y—  x)e   j  cosh  (kx)  e 

,  /_1  , • ,%      2.     cosh  (kL)  a*-  -X.) 
cosb(Y  -1  L)(X   -  aK  )  j 

2  2  (A-15) 

co  nnx      ALZ 

cos  —   e 


I 


2  2         2  2  2 

n  n  a  a.  x  wn  n  a  a.   2^2L_.  ..  .nn 


n=l,3,5...  ( — ,   +K«)( o  +  a*  M— )  sin(— ) 


Considering  the  second  term  in  equation  (A-10),  the  inversion 

as  obtained  through  use  of  equation  (A-ll)  is  given  in  equation  (A-16) 

s  s  ■  2 

n   q   a  A  -         „   q   a  A     \  t       a*  t 


j>-\       V       Hoou  "j         1       _V   Moou  "j  r  e  r     e^l_l_WA  .,, 
<\     L    k      (S-X  )(S-aK*)  ""  L     k   '^(X  -qk2)   (aK2-X  )  } 
j=l  J         j=l  j  J 


In  a  similar  manner  the  first  term  of  equation  (A-10),  the 
inversion,  again  using  equation  (A-ll),  is  given  in  equation  (A-17). 

z'1  j£  CjL cosh/F"  |=  p  Bi coah(y^ie Bit 


1  *  (S-fi.)cosh  /-  L     i=l      coshrV—  LI 
^i     \/  a  v   a  ' 


i,2    ' 

n^x       4L 
B   cos   2L   e 


2  2 

n  it  a  .  (A-17) 


I    I 


,,v..(.2LLfl  .  p  )  _L_  (Sln  BE) 

,  z     ■  i   niTd       2 

The  temperature  distribution  in  the  fuel  plate,  0f(x,t)   is  obtained 
*\>y   Substituting  equation  (A-15)  in  equation  (A-12)  and  adding  equations 
^A-12),  (A-16)  and  (A-17).   The  result  is 
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ef(x,t)=£  Bi«"V^»   >Plt   _y 


2    2 

-n  tt  a  t 

^^  oo  cos     -tt—  4L 

2L        e 


i^l      cosh    /'   .   .  -i^o    c      ,    L  nns 

'—  L  n-l,3,5,—  (— — )  sin(— ) 

a  anrr  l 


B,  ^     q      a  A      cosh  k  L                               . 

\  )     i +  V  -22 I y 

L_.         2    2  L             2    2                       2    2                .        J 

i-1     ILJLfi+p  J=lk(njEJl+x)(lLnJS+alc2) 

4L      1  4L      J    4L 


-   V  Ai  qoo  a  e  j    I  cosh(KL)  cosh(  /  J.  x) 

/_,  -1 J  cosh  (kx) j^=~ > 

j=l  k(aK  -  \   )        \  cosh(  J_i   L)     J 


(A-18) 


a 


As  seen  by  comparing  (A-l)  and  (A-2)  the  differential  equations 
to  be  solved  in  the  fuel  and  moderator  are  almost  the  same,  the  only 
difference  being  in  the  heat  generation  term.   Thus,  the  total  differential 
equation  for  the  moderator  in  the  Laplace  transform  domain  after  having 
applied  the  zero  Initial  temperature  condition  is  given  in  equation  (A-19) 

s 

2  S  V   F   A 

v   *ra(*'9)  -;♦.<«.»  =  - L  rrrfT  (A"l9) 

i=i  k  (S"V 

The  homogeneous  solutions  are  the  same  as  before  and  particular 
solutions  are  easily  obtained,  as  before,  from  the  method  of  undetermined 
coefficients.   Therefore,  the  solutions  to  equation  (A-19)  in  the 

transform  domain  are 

s 

F  a  A4 

(A- 20) 


$m(x,S)  =  C   cosh  /^\   +  D  sinh  /|  X,  +  £    F  a  Aj 


k  S(S-\  ) 
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The  transform  boundary  conditions  for  the  moderator  are  same  as 
those  for  the  fuel  if  the  origin  in  the  moderator  is  taken  at  the 
outside  of  the  unit  cell,  i.e., 

dC  (x,  ,  S) 
'Pi 


i=l 


and 


d  x, 


=  0  .    (A- 21) 


x,=0 


Thus  D  equals  zero  and  A  takes  the  same  form  as  for  the  solution 
in  the  fuel.   The  complete  solution  in  the  transform  domain  is  given 
in  equation  (A-22). 


4  (x,S)  =  ) 
m        L> 

1=1 


B,  cosh  /—  x. 

1 V  a  ' 

(S-p^cosh  y^L,    j=l 


E 


a  F  A 


k  S  (S-X  ) 


cosh  /—  x, 

[  1  -  ^-   1  . 


:°Sh  n/FL' 


(A-22) 


The  inversion  of  this  solution  is  easily  accomplished  by  the 
same  method  used  for  the  fuel  region.   The  solution  for  the  temperature 
is 


P  B.    coslW-1  x.fcPi*     °^ 

e  (x,t)  =  )     -^ 

m  L> 


nnx, 


2    2 
■n  tt  a 


cos    2L,         e      AL, 
1=1  cosh(v/^  L.)  n=T,3,5/"(    ^  >  8l"  ^~f~^ 


I 


il  -k — + 1 

vi-l        n  tt  a   ,    D  ,_. 


a  F   A, 

L 


2    2  2    2 


< 


+  Pj.        J=: 


n  tt  a       n  n  a 


} 


^ 


al; 


2     +ki 


8  PA  "KXt 

aFA  e   j 


k\ 


{ 


cosh 


1    - 


jT__x, 


cosh 


! 


L, 


(A-23) 
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This  solution  could  have  been  obtained  from  the  solution  in  the 

fuel  by  setting  k  =  0  and  q  =  F. 

"~        _po  "~ 

The  equivalent  solution  in  cylindrical  geometry  (r  dependence  only) 
for  the  exponential  boundary  condition  is  obtained  in  the  same  general 
manner,  however,  several  important  differences  do  occur.   The  equations 
governing  the  temperature  in  the  fuel  and  moderator  are 

r1'  qon  Xn  (*r)  A1  e  f  1    30f(r,t) 

V?e£(r,t)+[   °°  B      1 =1   -L-         CA-24) 

1=1 

and  * 

V.     F  A.  e   j     .90  (x,t)  ._. 

V2e(r,t)+V     3 =  1     -V  (A-25) 

m       / .        .        a     at 

j=l 
Considering  the  fuel  region  first,  transforming  with  respect  to 
time  and  applying  the  zero  initial  condition  yields 

^7  C(r,S)  -  -N(r,S)  =  -  )   — L- -2 (A-26) 

a  j=l   k  (S  "  V 

The  spatial  ooerator  for  this  case  is 

V2  =  -4-  +  !^  (A-27) 

a  2     r  or 

or 

The  oarticular  solution,  again  by  the  method  of  undetermined  coefficients, 

is 

q   a  A.  I   (Kr) 

<t  (r,S)  =  — ]   °  (A-28) 

P        k(S-\  )  (S-a-K  ) 

The  homogeneous  solution,  noting  the  fact  that  <£(o,S)  if  finite  is 

<£.  (r,s)  =  A  I   /?  r  (A-29) 

h  o  n/  a 
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The  final  constant  A  is  evaluated  by  use  of  the  surface  temperature 
boundary  condition  and  the  solution  in  the  transform  domain  is 


P    B,     I   /-  r     S„      q   a  A.      / 
r,s)  =  V  -* "k—  *  Y    -= L-T  (l 


R)  /2, 

n  Kin 


(kD 


o  v  a  o 


(A-30) 


i    /2,         / 


ova  o 


Inverting  this  expression  is  entirely  analogous  to  the  inversion 
of  the  equivalent  expression,  equation  (A-9)  of  this  appendix.   The 
analogy  carries  even  to  the  point  that  in  determining  the  zeros  of 
I   /—  R  as  in  determining  the  zeros  of  Cosh  /—  L  imaginary  values  of 
S  yield  an  infinite  set  of  zeros.   In  this  case  the  zeros  are 

I       /-  R     =   I       /—  R     =  I   1     /£  R     =  J       Mr     =0  (A-31) 

o\/ao         ova       o         ovao         ovao 

therefore, 

—  R  =  u  ,  where  w  =  2.4048,  5.5201,  etc.  (the  zeros  of  J(o  )=0) 
a  o.    n        n  n 


and 


^R2 

a     o 

2 

=  u 
n 

K  ■■ 

2 

W  a 
n 

R2 
o 

S 
n 

2 

u     a 

n 

R2 
o 

(A-32) 


The  details  of  the  remainder  of  the  Inversion  are  not  included 
The  result  is 


.         R    B.   I  (  /PI  rJePi1  „ 


2 

01   r  (,i  a 

z    n_  N       -     n 

J      (F~ }    e     F"' 
o 


1=1  x    t/Li     R   )  n=  1,2,3/ 


R2 


a  2oi  a  n 


B.  '  i  q       a  A.    I     (KR  ) 

i  '   \  oo  100 


HE   -r-4 — +'l       f    ]  °2  ° — } 

-T-  +.|3,  k   (-2—  +  V(~T  +   ^  } 


i        A.    q        a  eXjt      _  [I    (KR   )    I       /S.  rl 

-E  J  r  ,  o  { [I° (yr)] '  ~ — ^-z- 

»  a 


where  the  1*1  's  are  the  roots  of  the  equation,  JQ(X)  =  0. 

The  solution  for  the  cylindrical  geometry  in  the  moderator  is 
complicated  only  by  the  fact  that  the  symmetry  cond  i  t  io  n  cannot  be 
located  at  the  coordinate  r  =  0.   Thus  both  terms  I   /—  r  and 


o  *j  a 
K   /—  r  of  the  homogeneous  solution  must  be  retained.   The  details  of 

this  solution  are  not  included.   The  result  is 

e.«..o  -  f  b,  r  \/|F  »■>  W?  V  ;  S  ^S  r)         ,    eslt 


1   a 
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2     /V     [KJ^R)    Io(   7SV   I,  /^R,Ko(  /£,,].  V 


I 


+  >  f  a r__a L_a 

n=l 


X 


HvTS  rju/S  i9-K.ryS  «V?ji  ^io(7^  ^(VS  VoJ^  ^yi^*,] 

Tq  a  *  a  *  a  "a  r   a  fa 


P  B  s  S  a   F    A      eX1fc 

TV  i  V         a  F   Aj  |  +y  Q        Ai    6   3  (A-34) 

\k        ^      'A      -n<PnV  A 


X 


K,<V  -1  R,)I   («/-*■  r)+   I.t/-1  R,\K  (V-1  r) 
'7a       •      o  J  a  *  /   a       "  o  r   a 


Mv^MVy^'W^  v.'y'j1 


where  the  p   s  are  the  roots  of  the  equation, 

KY/^Y)I(    /^RJtl/  /^R,1K  C  /^Rj=    0 
,  lya     ''  o1  v  a  vy  a     '   o1  \/  a 

The  next  solution  considered  again  uses  an  exponential  fit  for  the 
heat  generation  rate,  however,the  boundary  condition  was  that  the  first 
derivative  with  respect  to  x  evaluated  at  the  outside  of  the,  plate  or 
effectively  the  heat  flow  out  of  the  plate  could  be  expressed  as  a  sum  of 
exponentials.   The  general  solution  with  the  exception  of  the  evaluation 
of  the  final  constant  A  is  exactly  the  same  as  the  first  derivation  in 
this  appendix.  Including  the  one  undetermined  coefficient  the  solution  is 


rs  r1       qoo  a  Ai  cosh  yX 

<E,(x,3)    =    A  cosh    /-  x     +  )  ] —      .  (A-35) 

...j        k(s-A.,)    (S-ctK  ) 
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Evaluation  of  A  through  use  of  the  boundary  condition  leads  to 


^  Bi  c°sh<y^x)      v  q°oa  ai 

$f  (x,  S)  =  )  -^ *-S +  ) ] — - 

i=l  (S-b.)  /-  sinh(  /—  L)  1=1       j 

'  i  n/  a     V  a   '  J        J 


,.„«i,  v     cosh(KL 
cosh  kx  -  


(A-36) 
)coshw  a  x)1 


oshw  a  L)1 


The  details  of  this  inversion  are  not  included.   The  result  is 


2   2 

n   tt   x 


P     /    Bt   cosh(y/-x)ePi  .B^  cos(-2TT)e 

M    /Pi  . .  /Pi .   "  PiLi"  v 


1=1*        n/^  sinhV^L)  ^    n-l"  a  ( 

a  K     a        } 


y    J    V  Bi  V  qoo   a   A1    ^  Slnh  KL  I  (A-37) 

X    {    L  2    2  +     L  2    2  2    2  _        J 

4-1      n  n  a       d  4-1  wn  ^  a.   ,    ,    x  ,n  tt  a   ,         2.      J 

i-1     — ^ —  +  8.       J-i  k( r-  +  \   )( —  +  a*  ) 

4L       1  4L      J   AL 


X  A? 

£'   q   a  A,  e   j  v  sinh(XL)cosK  v    x) 

\    oo 1 r   ,      a    -, 

>    ^ [cosh  Kx  -  7==* = ]   -  q   a  A 


k(aK  -  M 


J-*-   slnh  V  -J-  D 
»  a      Y  a 


sinh  kL 


4=1   "xv**     -4 

k  \  KL 


The  time  dependent  temperature  distribution  in  the  moderator  for 
the  comparable  boundary  conditions  is  obtained  by  setting  %   equal  to 
zero  and  q   equal  to  F  in  equation  (A-37).   The  result  is 


P'  /  B  cosh(y  —  x)e  Pifc      B^ 
Mx,t)  =  I  — — -   °   mj.         -p-T 

i=l  »    J—  sinhrJ—  U 


cy-M 
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n   T1    a 
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oo                           nTTx,  P  r 

■I         °-^^ (I  -TT±—} 

n=l                     I  -  J  cos   nn  i-1               —  +   R 

a  4L2  * 


8' 


(A-38) 


a  F   A,  \.t        ,  % 
L        (  e     j      -    1  ) 


j=i    k(V 


The  final  solution  considered  again  used  a  sum  of  exponentials  to 

represent  the  time  dependence  of  the  heat  generation  rate,  however,  the 

surface  temperature  boundary  condition  was  approximated  by  an  even 
P 

Fourier  series,  )        B  cos  R.t.   With  the  exception  of  the  steady  state 

1=1 
term  resulting  from  the  surface  temperature  boundary  condition  the 

derivation  follows  exactly  the  first  derivation  of  this  appendix.   The 

steady  state  term  is  handled  most  easily  in  a  slightly  different  manner 

as  shown  below.   The  general  solution  in  the  transform  domain  is 


$A 


P      '  n    o  t_      /S  o  q        a   A . 

B.S  cosh/   /-  x-v  *  Moo  j 


^ 


,_.    rc2,R2.  /?*  £,         k(S-\.)    (S-qk2) 

i-1    (S  +p.)    cosh(\/  a  L)        J-1  j 

/$"*  (A-39) 

a     x) 

[COSh    KX    -    "  -rp  ] 

cosh(y  a     L  ) 


The  first  term  presents  the  only  change  and  there  only  for  the 
poles  at  S  =  +  jR .   The  terms  generated  from  the  inversion  integral  by 
these  two  poles  are 


|       B   j  B   cosh(>/—  x) 
P,      = i 7?§= e  JPiC  (A-40) 

I         2  j  R    oshCy/— i  L) 
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and 


B     (-jp   )   coshp/—  x) 

=     r~  Mr>-\ e        *  i 


2    (-jp 


1)    cosh(Y  ~ —  L  ) 


(A-41) 


S=- 


iPi 


These  terms  are  most  easily  handled  by  recognizing  the  fact  that  the 
sum  p.  +  P9  is  the  sum  of  a  function  and  its  conjugate.   That  is 


p.  +  p  =  f(z)  +f(z)=  f(z)  +  f(z)    since  z  is  a  pure  imaginary  number, 


Therefore 


Pt  +  P2  =   2  Re  |  f(z)   }  =    2  Re|p1  | 


Bt  coshf^-^- x)      .  . 

7f7~  e       J 

2  cosh(/s/ — £  L) 


(A-42) 
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i      e 


cosh^  — -  xj 


cosh 
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cosh 


cosh 


B     Zj|    (x)      cos   [pit  +  rp1(x)] 


where 


Z    (x)   = 


cosh 


^i 


cosh  *J  — -  L 

a 


cosh(y  ~  x)cos(y-^  x)+    j    sirxy— x)  sinh(y  —  3$ 


[tj  T-  Li+    1    sinrv7  —  L)  sinhrJ  -r- 


/  ~  L)cos(y  —  L)+    j    sin(y  ^  L)  sinh(y  —  I 
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"IT 

2a 


cos 


2^X>] 
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(A-43) 
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The  resultant  time  dependent  temperature  distributions  are 


2   2 

n  tt  a 


00  nnx  4L/ 

cosi>L    e 


e£(«.«  =  ^  b,  zj*  c08  (p  t  +  n>  -  I 

._,  „_,    ,    c   ...(L   /nrra)      sin  — 

i-1  n-1,3,5,  2 


nir 


'{ 


?    Bi 


2    2 

,n_n_cu 

2 

4L 


i=l      ( 


4  4   2 
iL^t-^12) 


Jt,       q       a   A     cosh  kL 


16L 


2    2  2    2 

-1      k^-V^  aK2)(ILJlJX  +  X.) 
4L  4l/  J 


(A-45) 


•     )  r— ■ -1    COSh 

^       k(aK     -  \  )  l 


KX 


cosh  (kL  )  cosh(y  — L  x) 


cosh(V  — L     L  \ 

v'   a  ' 


87 


and 
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in  the  fuel  and  moderator,  respectively 
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APPENDIX  B 


Description  and  Explanation  of  the  IBM-650 
Computer  Program  Used  for  Fitting  Empirically 
Experimental  Data  with  the  Sum  of  Several 
Terms  of  Exponential  Form 


The  computer  code  was  written  to  fit  an  analytical  function  of 
the  form  of  the  sum  of  exponentials  to  the  experimentally  determined 
power  traces  during  a  transient  burst.   The  program  was  written  in 
SOAP  II  and  floating  point  form.   The  object  program  is  listed  and 
the  logic  diagram  is  shown  in  this  appendix. 

The  criteria  that  the  machine  inspected  was  that  the  sum  of  the 
squares  of  the  residuals  between  the  experimental  data  and  the  calculated 
values  should  be  made  as  small  as  possible.   Each  of  the  fitting  para- 
meters was  varied  in  turn  by  a  specified  increment,  holding  all  other 
parameters  constant,   until  such  a  time  that  a  specified  increment 
could  make  no  further  reduction  in  the  sum  of  the  squares  of  the  residuals 
This  parameter  was  then  stored  as  the  best  available  estimate  of  the 
particular  empirical  parameter.   When  none  of  the  parameters  could  be 
varied  by  the  specified  increment  to  give  a  smaller  sum  of  the  squares 
of  the  residuals,  the  increments  were  refined  and  the  trial  and  error 
process  was  repeated  with  the  refined  increments.   This  procedure 
continued  until  the  increments  were  less  than  a  specified  precision. 

The  data  were  fit  empirically  with  a  function  of  the  f~rm 
s 

V±  =   V  A.  e  Xjtl  s  =  10  (B-l) 

j=l 

As  stated  above,  the  best  fit  criterion  was  that 

c  s 

:  j  i   -  <t,  )  '  (B-2) 


-I  ^  !(vjl  -v 


1=1  *    j=i 
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be  a  minimum. 

The  program  required,  in  addition  to  the  experimental  data  and 
their  respective  times,  initial  estimates  for  a  specific  number  of 
parameters.   The  program  had  a  capacity  for  up  to  20  data  points  and 
10  terms  in  the  summation  of  equation  (B-J.)  .   These  input  data  were 
read  into  the  machine  along  with  the  program  deck  on  one-word  load 
cards.   Each  one.-word  load  card  contained  a  particular  constant  or  an 
initial  value  and  its  specified  storage  location.   Table  B-l  lists 
the  various  input  data  needed  for  this  program. 


Table  B-l.   Input  data  required  for  use  of  the 

IBM-650  program  which  fit  empirical- 
ly experimental  data  with  several 
terms  of  Exponential  form. 


Symbol 


Explanation 


Storage  Location 


ZERD 

FPONE 

FOUR 

TEN 

HNDRD 

ONE 

EIGHT 

INDXB 

INDXA 

ONEHD 

DELBj 

DELRj 

BINITj 

RINITi 

Zi 

FLUXi 
W. 


0.00 

1.00 

4.00 

10.00 

100.00 

Index  Number  1  (0000000001) 

Index  Number  8  (0000000008) 

Number  of  Exponential  Terms  (OOOOOOOOxx) 

Number  of  Data  Points  (OOOOOOOOxx) 

Precision 

Initial  increments  of  the  Amplitude,  B, 

Initial  increments  of  R< 

Initial  estimate  of  Bj 

Initial  estimate  of  Rj 

Time  of  ith  data  point 

Experimental  data  at  ith  point 

Weighting  function  at  ith  point 


J 


0073 
0129 
0185 
0158 
0090 
0392 
0024 
0412 
0062 
0168 

(1800  +  i 
(1810  +  j 
(1300  +  j 
(1310  +  j 
(1200  +  i 
(1200  +  i 
(1750  +  i 
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The  machine  yielded  an  answer  card  having  a  capacity  of  8  words, 
a  word  being  ten  digit  numbers  and  a  sign.   For  the  first  answer  the 
machine  punched  out  the  initial  estimates  of  the  fitting  parameters  on 
as  many  cards  as  was  necessary  to  accommodate  them.   A  and  \  for  the 
first  term  were  stored  in  word  locations  1  and  2,  respectively.   After 
a  card  was  filled  to  its  8  word  capacity,  it  was  punched  and  a  new 
card  began  to  fill.   This  procedure  continued  until  all  of  the  fitting 
parameters  had  been  punched  out.   Then  a  separate  card  was  punched 
giving  in  the  word  8  location,  the  value  of  the  sum  of  the  weighted 
of  the  residuals  between  the  experimental  data  and  the  calculated  values, 
The  machine  then  punched  out  values  of  the  time,  residual  and  correct 
values  at  the  last  data  point  in  word  locations  1,  2  and  3,  respectively 
and  the  next  to  last  data  point  in  word  locations  5,  6  and  7,  respec- 
tively.  The  same  information  for  two  previous  data  points  was  punched 
out  on  a  second  card  in  the  same  format  and  this  procedure  was  continued 
until  the  position  residual  and  correct  value  was  punched  out  for  each 
data  point.   Subsequent  improvements  in  the  parameters  and  the  weighted 
sum  of  squares  of  the  residuals  were  printed  out  after  each  cycle  of 
trying  to  vary  each  parameter.   The  positions,  residuals  and  correct 
values  at  each  data  point  were  obtained  at  this  time  if  the  console 
instruction  was  negative.   When  the  fitting  parameters  could  not  be 
further  improved  with  the  most  refined  increment  specified,  the  punching 
of  the  best  fit  parameters,  the  sum  of  the  weighted  square  of  the  resi- 
duals, the  position,  residual  and  correct  values  took  place  according  to 
the  procedure  described  above. 
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OBJECT     PROGRAM-APPENDIX     B 


Q2 


REtIT 

8  TUUF 
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APPENDIX  C 


Description  and  Explanation  of  the  IBM-650 
Computer  Program  Used  for  Fitting  Empirically 
Experimental  Data  with  an  Even  Fourier  Series 


The  computer  code  was  written  to  fit  an  even  Fourier  series  to 
the  experimentally  determined  surface  temperature  traces  during  a 
transient  burst.,  The  program  was  written  in  SOAP  II  and  floating 
point  form.   The  object  program  is  listed  and  the  logic  diagram  is 
shown  in  this  appendix. 

The  data  were  fit  empirically  by  a- finite  number  of  terms  of  the 

even  trigonometric  series. 

P 
6<t)      =  V       B,      cos  lllii 


i  a 

1=1 


(C-l) 


where 


B     =   2/ a  ja  y(t)   dt 


and  B,   =   1/a  Ja  y(t)    cos     ^^     dt    . 

i  o  a 

The  integrations  were  carried  out  numerically  by  means  of  Simpson's 
rule  thus  requiring  an  odd  number  of  data  points. 

The  program  input  consisted  of  the  experimental  data  and  their 
respective  times,  the  period,  the  time  increment  between  data  points 
and  a  specification  of  the  number  of  terms.   These  data  were  read  into 
the  machine  on  one-word  load  cards.   Each  one-word  load  card  contained 
a  particular  constant  or  piece  of  data  and  its  specific  storage  location. 
Table  C-l  lists  the  input  data  needed  for  this  program, Storage  locations 
limit  the  product  of  the  number  of  terms  and  the  number  of  data  points  to 
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less  than  550.  The  number   of  terms  and  the  number  of  data  points  are 
each  limited  to  less  than  50. 


Table  C-l.   Input  data  required  for  use  of  the 

IBM-65Q  program  which  fits  empirically- 
experimental  data  with  a  finite  number 
of  terms  of  an  even  Fourier  series. 


Symbol 


Expanation 


Storage  Location 


ZERO  0.00 

FP0ME  1.00 

FFTW0  2 .  00 

FFTPvE  3.00 

FPF0R  A. 00 

PI  3.14159 

INDX2  Index  Number  2  (0000000002) 

INDXJ  Number  of  Terms  (OOOOOOOOxx) 

INDXK  •  Number  of  Data  Points  (OOOOOOOOxx) 

A  Period'of  Cosine  Terms 

H  Time  Increment  Between  Data  Points 

Ti  Time  of  ith  data  point 

Yj_  Experimental  data  at  ith  point 


0083 

0034 

0108 

0207 

0008 

0261 

0160 

1001 

1002 

1003 

1004 

(1100  +  i) 

(1150  +  i) 


The  machine  punched  out  a  card  having  an  eight  word  capacity,  each 
work  consisting  of  10  digits  and  a  sign.   The  first  output  consisted  of 
the  time,  the  calculated  value  and  the  residual  between  the  calculated 
values  and  the  experimental  data  for  the  last  data  point  in  word  locations 
1,  2  and  3,  respectively.   The  same  information  for  the  next  to  last  data 
point  was  punched  out  in  columns  5,  6  and  7  of  the  same  card.   The  same 
information  for  the  two  previous  data  points  was  punched  out  on  the  next 
card.   The  above  procedure  was  continued  until  the  time,  residuals  and 
calculated  values  were  punched  for  each  data  point.   B.  was  then  stored 
in  word  location  1,  B_  in  word  location  two,  etc.  until  all  of  the  B's 
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had  been  stored  and  punched.   If  more  than  8  B' g  were  calculated  the 
fir9t  8  were  stored  and  a  card  punched.   Additional  B' s  were  punched  in 
succeeding  cards  with  the  lower  number  B' s  starting  on  the  left  of  each 
card.   Finally  the  summation  of  the  square  of  the  residual  at  each  data 
point  was  punched  out  in  word  location  8  of  a  final  card. 
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OBJECT  PROGRAM-APPENDIX    C 


100 


E  0  OC  A 

N  E  C  •  T 


R  E  OH  C 


C  08  I  0 


E  0  OSR 

HECK 

R  f  nu  u 

8  I  NET 


NCG9I 


E  nuf  f 
line 

T  ■  OP  I 

ONEPI 
F  P  ONE 
COS 


I  NTCT 


L  OP  33 


0  0  1  1  1 


E  0  OC  L 


zmo 

START 


SLA 

8  Y  N 
BIB 
ITN 
8  Y  N 
8  Y  N 
8  Y  N 
8  Y  N 
8  Y  N 
8  Y  N 
8  Y  N 
8  Y  N 
8  Y  N 
8  Y  N 
8  Y  N 
8  TO 
8  U  I 
F  »  0 
BU  I 
FSB 
FSB 
BU  I 
F  A  0 
S  TU 
RSI) 
S  T  U 
8  T  U 
STL 
8  TO 
BM  I 

f  a  o 

BM  I 
FSB 
FSB 
BU  I 
F  A  0 
8  TU 
R8U 
8  TU 
S  T  U 
LOO 
8  TO 
R  A  U 

r  ao 

8  T  U 
F  A  0 
8  T  U 
R  S  U 
FMP 
FM  P 
F  0  V 
F  0  Y 
S  T  U 
RAM 
STL 
RAM 
R  A  U 
F  0  Y 
FSB 
BU  I 
R  A  U 
F  A  0 
8  T  U 
R  A  U 
10 
6  2 
31 
10 
S  T  0 
R  A  U 
FMP 
LOO 
S  T  U 
S  TO 
R  A  B 
R  A  U 
S  TU 
A  X  8 
R  A  U 
FMP 
F  A  0 
S  T  U 
A  «  B 
R  A  U 
FMP 
FAD 
8  T  U 
R  A  U 
SUP 
A  U  P 
N  I  U 
A  X  B 
R  A  U 
FMP 
F  A  0 
S  T  U 
A  X  B 
R  A  U 

r  a  o 

FMP 
F  0  V 
8  T  0 
LOO 
8  TO 
8  T  0 
8T0 
8  T  0 
8  T  0 
S  TO 
8  TO 
S  T  0 
00 
R  A  U 
UP  Y 
STL 


10  0  0 
ERROR 
09  7  7 
T 
Y 
8 

8  T  A 
P  H  I 
FLUX 
C  0  S  B  T 
START 
I  N  0  X  J 
'  I  N  0  X  K 


EXIT 

N  E  C  A  T 

T  W  0  P  I 

N  E  G  A  T 

0  NE  P  I 

T  »  0  P  I 

ON  E  P  I 
T  H  E  T  A 
F  P  0  N  E 
T  E  RM  M 
F  U  N  K  T 
E  N  N 
EXIT 
N  E  C  A  Y 
T  *  0  P  I 
N  E  C  A  Y 
0  N  E  P  I 
TIOPI 

0  N  E  P  I 
T  H  E  T  A 

B  00  3 
T  E  RM  M 
F  U  N  K  T 
F  P  0  N  E 
E  N  N 
E  N  N 
F  P  ON  E 
N  P  0  N  E 
F  P  0  N  E 
E  N  N 
T  E  RU  U 
T  H  E  T  A 
T  H  £  T  A 
N  P  0  N  E 
E  N  N 
T  E  R  MM 
F  U  N  K  T 
F  M  A  G 
T  E  RUM 
8  00  3 
F  U  A  G 
SIZES 
E  N  U  F  F 
F  U  N  K  T 
T  E  R  U  M 
F  U  N  K  T 
F  U  N  K  T 
00  0  0 
8  3  18 
4  15  9 
0  0  00 
E  X  I  TC 
B  T  A 
T 

C  0  S  B  T 

E  G  S  I  T 

00  0  1 
P  H  I 

A  C  C  U  M 

00  0  1 
P  H  I 

F  P  F  0  R 

A  C  C  UM 

A  C  C  UU 

00  0  1 
P  H  I 
FPtIO 

A  C  C  U  U 

A  C  C  UM 

8  006 

1  N  0  X  K 
I  N  0  X  3 
L  OP  33 

0  0  0  1 
P  H  I 
F  P  F  OR 
A  C  C  UU 
A  C  C  UM 

0  0  0  1 
P  H  I 
A  C  C  UM 
H 

F  P  TR  E 
Z  Z  Z  1 
Z  Z  Z  10 

09  7  7 
09  7  8 
097  9 
09B  0 
09  B  1 
09  8  3 
0  9B  3 
09  6  4 
00  00 

I  N  0  X  K 
I  N  0  X  J 
I  N  0  X  C 


1»99 

10  SO 
098  5 
1100 

11  SO 

1  aoo 
laso 

1300 
13S0 

1400 
1000 
1001 
1003 
10  03 
1004 

«moc 


cos  i  o 

R  E  0  U  C 
CO  t  I  0 


NCO  ST 

RCOUD 

S  I  HE  T 

RCOUO 
8  I  NE  T 


NEO  ST 


NEC  ST 
EXIT 

004  3 

5  35  1 
37  5  1 
00  5  1 


E  00  CH 
E  X  I  TC 


L  OP  J  3 


G  0  1  1  1 


E  GS  I  T 


Z  Z  Z  1 

00  00 


1  0000 

3  00  0  0 

3  0  0  00 

Tl   SMALL         4  0000 

Yl   WITH   Tl      5  0000 

6  00  0  0 

7  00  0  0 
P  000  0 
9  0  0  0  0 

10  0  0  0  0 

11  0000 
NO  TERMS  12  0000 
NO   OATA   PT    13  0000 

14  0  0  0  0 

15  0000 

16  00  0  0 

17  0006 
IB  00  0  9 

19  00  39 

20  00  4  3 
3  1  0  0  10 
3  3  0  0  8  9 
3  3  00  4  3 
3  4  00  3  3 

25  00  3  1 

26  0  139 

27  00  4  7 

28  000  5 
3  9  00  5  0 

30  0056 

31  0  109 
3  2  0169 

3  3  009  3 

34  00  6  0 

35  0  3  39 

36  00  9  3 

37  0073 

38  00  8  1 

39  0  3  8  9 

40  009  7 

41  005  5 

4  3  00  3  7 
4  3  00  6  3 
4  4  00  13 
4  5  00  11 
4  6  00  19 
4  7  00  6  1 
4  8  0  113 

4  9  00  4  9 
50  00  7  8 

5  1  0  13  8 
5  2  00  6  6 
5  3  0  159 

54  014  7 

55  000  7 

56  00  14 

57  00  99 

58  0  0  5  7 

5  9  0  16  1 
60  00  4  1 

6  1  00  4  5 
6  2  0  10  7 

63  00  3  1 

64  009  4 
«5  006  4 

66  00  13 

67  00  4  6 

6  6  00  3  4 
«  9  0  10  0 

70  0  10  6 

71  0  10  5 
SlIBLOSUP       73  0150 

7  3  0  10  3 

74  0  300 

75  0156 

76  0  16  3 

77  015  5 

78  006  3 

7  9  00  6  9 
INPUT            80  0205 

8  1  0  0  5  8 
8  2  0  0  8  7 

83  0  113 

84  0  119 
INPUT           85  0255 

86  01  S  8 

87  013  7 

88  016  3 

89  007  1 
INPUT        90  0157 

91  00  15 

93  00  30 

93  00  36 

94  0  30  5 

95  0  20  8 

96  0  18  7 

97  0313 

98  0  16  9 

99  0  3  5  5 
INPUT  100  0237 
INPUT  101  0004 
EXIT   INSTR   103  0350 

10  3  0  20  6 

10  4  0  3  13 

10  5  00  3  0 

10  6  0  13  1 

10  7  00  3  3 

10  8  0  0  3  3 

10  9  0  0  8  4 

110  00  3  5 

111  0036 
113  0  3  0  9 

113  100  0 

114  0  3  5  7 

115  00  3  3 


00  00  00  0000 

00  0000  0000 

00  0000  0000 

00  0000  oopo 

00  0000  oooo 

00  0000  0000 

00  0000  0000 

00  0000  0000 

00  00  00  0000 

00  00  00  0000 

00  0000  0000 

00  0000  0000 

00  0000  0000 

00  00  00  0000 

00  0000  0000 

34  0003  0006 

46  0009  0010 

33  00  13  0039 

46  00  09  0043 

33  00  46  0033 

33  00  12  0089 

46  0043  0010 

33  00  46  0033 

21  00  38  0031 

61  00  34  0139 

31  0044  0047 

31  0003  0005 

30  0059  006? 
24  0003  0056 
46  0109  0060 

32  00  13  0189 
46  0109  0013 

33  00  46  0073 
33  00  12  0239 
46  00  93  0060 

33  0046  0073 
21  0020  0001 
61  8003  02H9 
21  00  44  0097 
21  00  03  0055 
69  00  34  0037 

34  0059  0063 

60  00  59  0013 

33  00  34  0011 

31  00  16  0019 

32  00  34  0061 
21  00  59  0113 

61  00  44  0049 
39  00  36  0078 
39  003b  013B 

34  0016  0066 
34  00  59  0159 

31  00  44  0147 
67  00  03  0007 

30  0111  0014 
67  00  44  0099 
60  8003  0057 
34  0111  0161 

33  00  64  0041 
46  0094  0045 
60  00  03  0107 

32  00  44  0031 

31  00  03  0063 
60  00  03  0003 
10  0000  0043 
63  8316  5351 
31  4159  37  5  1 

10  0000  0051 
3  4  '00  5  3  0  10  6 
60  3350  0105 
39  5100  01  SO 
69  0103  0000 
31  7400  0053 
24  0153  0156 
63  0001  0162 
60  5300  0155 
31  0110  0063 
53  00  01  0069 
60  5300  0305 
39  0008  0058 

33  0110  00B7 
31  0110  0113 
53  00  01  0119 
60  5300  0355 
39  0106  01 5 B 
33  0110  0137 
31  0110  0163 
60  8006  0071 

11  1003  0157 
10  0160  0015 
44  00  63  003C 
53  0001  0036 
60  5300  0305 
39  00  OB  030B 
33  0110  0187 
31  0110  0313 
53  0001  0169 
60  5300  03  5  5 

33  0110  0337 
39  1004  0004 

34  0307  015  3 
34  0303  0306 
69  0309  0313 
24  0977  0030 
24  OVT  U  0131 
24  0979  0032 
24  0980  0033 
24  09P1  0084 
24  0983  0035 
34  0983  0036 
34  0984  0303 
00  0000  0000 
60  1002  0357 
19  1001  0033 
30  00  37  0080 
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8 
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8  T  0 

09  8  7 

1 

315 

A  X  A 

00  04 

216 

8  X  B 

00  0  1 
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52 
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0  5  11 

0  511 

69 
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24 
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0  40  3 

53 

00  0  1 
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02  5  9 
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026  3 
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69 
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39 

0  10  8 
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59 
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34 
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21 
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0  60  3 

51 
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03  59 
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03  10 
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0313 
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0269 

82 

00  0  0 
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52 
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69 
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50 
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88 
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69 
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0  3  30 

5  8 
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6  9 
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0  2  5  4 

02  5  4 

80 
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69 
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0  3  8  9 
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0  16  5 

21 
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51 
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014  8 

0  14  8 
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52 
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APPENDIX  D 


Description  and  Explanation  of  the  IBM-650 
Computer  Program  Used  to  Calculate  Temperature 
Distributions. 


The  computer  program  was  written  to  calculate  the  temperature 
distribution  in  a  unit  cell  of  a  nuclear  reactor  system  given  the 
heat  generation  .rate  and  fuel  element  surface  temperature  as  a 
function  of  time.  The  temperature  rise  over  the  initial  temperature 


is  given  by 

n  tt  a 


0f(x,t)  =  V  B1  Z*  cos  (p£t  +   y±)    -  £ 


co    cos^— — )       e 


4L2 


t 


i=l  n=l,3,5/»-  (L  /nna)  (sin  j1) 


2  2 

..  1  f  Bi  (  ^~  )  ,   y     W  A1  cosh  KL t     (D_1} 

\L  4  4  2    9     Zj         2  2n      2  2 

1=1   ^4-  +  p2    j=l    k  (*-*=&*  )(njLS  +  aK2) 
16IT     1  4l/   J   4l/ 


^-,    q   a  A  e  j      coshfy  — *•  xjcosh  kx 
j-l  j  cosrUy  — L  LJ 


in  the  fuel  and  by  2  2 

_n  tt  g 

P       ,  °°     cos(M)  e   4L 

r-i           1  r-i      COSV      /   e 

9  (x,t)  =  )    B^  Z2  cos  (R.t  +  q>  )  +  )    — ^ 

m        Z_j    i   i  '  1    ri     Z_i     ,T  2  .    w    n^s 

1=1  n=l,3,5/^L  /nTTa)(slnr) 


t 
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2    2„ 

n         /  n   n   An 

4Z%r^-  -  .1    ^ ( 

i=i  ILj4j-  f  p,      j=i  k  c^^  +  x.) 

16IT  i  4l/  J 


(D-2) 


<-,         F  a  A.    e     j  r      cosh  J  — *■  x 

+  y   — j — -    {  — **—  - 1 

i=l  k   <-Xj)  cosh/^L 

in  the  moderator.   The  moderator  equation  is  obtained  from  the  fuel 
temperature  distribution  by  setting  q   equal  to  F  and  K  equal  to  zero 

The  equivalence  between  elements  of  the  algebraic  equations  and 
the  symbolic  logic  of  the  computer  program  is  shown  in  Table  D-l. 


Table  D-l.  Definition  of  symbolic  terms  of  the 
IBM-650  computer  program  for  calcu- 
lating temperature  distributions. 


Al 


.  =    B  Z  (x)  cos  p  t 


A3P.  =  q  a  A.  aX]t/   k  (ax2   -  X.) 


.„,,,,,,    cosh  kL  cosh  J — L  x 

A3SUM.  =  1   a         , 

j   7t==s -  cosh  KX 


cosh  si  —   L 


A3,  =    (A3P.)  (A3SUM.) 
CSHLL  =  cosh  J2±   L 


CSHLX.  =  cosh 
J 


COSLL  =  c 
J 


COSLX  =  cos 


Table  D-l  cont 


ARG1     =   nn   2L 
n 
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ARG2     =   n  n  a     4L 
n 

2   2 

n_n_cu 

A2STli  =  B      \h2 


4  4   2 

n  tt  a        .    r,    2, 


+ 


16L 


Pp 


A2ST1 


J 


a     a  A.    cosh  kL  /  k   ( 
co        j  '  , „  2 


2   2 

_I 

4L: 


2    2 


+   X    ) 


+       QK    ) 


2    2 

n  tt  a 


cos 


A2DDT     = 


2L 


4L' 


(-— -  )    sin     — 

nnct  2 


41/ 


P  s 

A2SUM  .  =    (   V     A2ST1     +  V     A2ST 
1=1    ,  j=l 


j   n 


fE.Rli  =  (A2DDT  )  (A2SUM  ) 
n         n        n 


Table  D-2.   Input  Data  Required  for  Use  of  the 
IBM-650  Computer  Program  Used  to 
Calculate  Temperature  Distributions. 


Symbol 


Explanation 


Storage  Location 


ZERO 

0.00 

ONE 

1.00 

TWO 

2.00 

PI 

3.14159 

FIFTY 

50.00 

CRIT 

0.0001 

ALPHA 

Thermal  Diffusivity 

0264 
0662 
0520 
0018 
0361 
0788 
0278 
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Table  D-2  cont. 


Symbol 


Explanation 


Storage  Location 


KAPPA 

Qoo 
KAY 

L 
INDXM 
NOLAM, 
AYEJj  ' 
LAMDAj 
AMMM-l 
BTAAi 


Reciprocal  of  Thermal  Neutron  Diffusion  0436 
Length  in  Fuel 

Normalization  Factor  for  Heat  Generation  0324 

Thermal  Conductivity  0581 

Half-Thickness  of  Region  0186 

No.  of  Terms,  Surface  Temperature  Fit  (OOOOOOOOxx)  0076 

No.  of  Terms,  Heat  Generation  Fit  (OOOOOOOOxx)  0456 

Amplitude  Parameter,  Heat  Generation  Fit,  Aj  (0200  +  j) 

Exponential  Parameter,  Heat  Generation  Fit,   j  (0220  +  j) 

Amplitude  Parameter,  Surface  Temperature  Fit,  B^  (0100  +  i) 

Period  Parameter  Surface  Temperature  Fit,  B^^  (0120  +  i) 


The  output  from  this  program  is  punched  out  on  one  card  having  an 
eight  word  capacity,  one  word  consisting  of  10  digits  and  a  sign.   The 
form  of  the  output  is  shown  in  Table  D-3. 


Table  D-3.   Output  Form  for  IBM-650  Computer 

Program  Used  to  Calculate  Temperature 
Distribution 


WORD  1  WORD  2 


WORD  3 


WORD  4 


WORD  5  WORD  6  WORD  7  WORD  8 


i=l 


CO 


6(x,t)   £ai4    I 


(Term) 


n=l,3,5 


I    A3j 

3=1 


33- 

o  •    • 
UtO  to 


id 

•  o 


to    <" 
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-1 

_l 

UJ  *~ 

CE 

Hx 

t 

(0 
Z 
< 

tO_J 
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< 

*  s» 

E 

<  o 

o 

f.  <"> 

H 

u. 

CO 

Z 

0.  X 

IL 

o 

K 

UJ* 
Vx 

ff> 

IK  ■ 

00  — 

K 

Z 

•  • 

O 

to  to 

u 

<M 

a. 

♦  ♦ 

\      °- 

ht 

J      o 

»-h 

_i 

U  UJ 

OTco 

El- 

S*-T 

.toi. 

Z  •* 

5oOTOT 

Sf~: 

gwto 

O  N(M 

d.  «< 

Til- 

■ s 

s'8 

t- 

T.X 

•  z 

K 

O   w 

to  — 

< 

<?£ 

c/i 

So" 

\z. 

•«B 

E 

« • — i 

»r 

*-*-- 

o  c 

I-^» 

5*  '" 

2 

a 

ito 

O 

a. 

oO  - 

111 

K   O 

o 

IL  U< 

CC 

-i 

(fljf 
to  ioio  — 

o  o  o  • 
t£  It.  U.  10 


«S\ 

Form  CSHLXj  8  CSHLL) 
Form    A3SUMj 
Form    A3j 
S*t  j-i-l 

fij 

to 

,g«*UI_N 

*°  .  cccc 

Z  Z  £  «j  ,, 
■    ■   «M 

*-«  Eg 

•  •loo 
wto  to  U.U. 


o 
111     (t 

Z       Ul 
O  Ul  M 

■  z  • 

<  °Q 

?zx 

•  «« 

to  to  to 


OBJfCT    PRO0RAM-APPEN0IX    U 


K1H 


HI  H 

o  a  no 

0H60 

HL  H 

1  Hno 

1  M  9  9 

SYN 

A  H  T  A  N 

1  U  00 

SYN 

ST  A  H  T 

1999 

8Y» 

A  Y  E  J 

o  a  oo 

SYN 

L  A  M  n  A 

0280 

SYN 

A  3 

0  2  4  0 

Bl« 

19  5  1 

19  60 

HI  H 

197  7 

19  B  4 

BLR 

0  1  OO 

0160 

SYN 

AMMO 

0100 

8  Y  N 

H  T  A  A 

0130 

SYN 

ALPHA 

087H 

SYN 

KAPPA 

0  4  3d 

SYN 

L 

01  Bfc 

SYN 

000 

0  324 

SYN 

KAY 

05H  1 

SYN 

1  N  0  X  M 

007  6 

SYN 

N  0  L  A  U 

0  4  5b 

SYN 

A  1 

0140 

PM  1  9H 

S  T  0 

•  BCD') 

R  A  U 

F  C  L 

FUP 

F  C  X 

SID 

D  F  N  N  1 

R  A  U 

F  SL 

FUP 

F  S  X 

F  A  0 

0  F  N  N  1 

8T1) 

0  N 

RIU 

F  SL 

FM  P 

F  C  X 

8  T  U 

N  11  M  M  a 

R  A  U 

F  C  L 

FUP 

F  S  X 

FSB 

M  0  M  U  A 

FO  V 

ON 

LOO 

utim 

811) 

IBEIO 

RAO 

ON 

RM  1 

A  B  C  02 

RAO 

THEIO 

F  A  0 

P  1 

3  T  0 

THETO 

A  B  C  02 

EOOE« 

S  T  0 

A  A  A  1 

3  T  0 

A  A  A  1  4 

RAO 

A  A  A  1  6 

F  A  M 

A  A  A  1  4 

3  T  0 

A  A  A  a 

LOO 

A  A  A  3 

3  TO 

A  A  A  4 

A  A  A  (A 

»  A  A  8 

RAO 

AAA.) 

r  g  b 

A  A  A  5 

»"  1 

A  A  A  6 

9T  U 

ill 

RAO 

A  A  A  4 

FM  P 

A  A  A  7 

3  T  0 

A  A  A  4 

A  A  A  (t 

A  A  A  0 

RAO 

A  A  A  2 

F  9  B 

A  A  A  3 

B  M  1 

A  A  A  28 

3  T  0 

A  A  A  2 

R  A  0 

A  A  A  4 

FUP 

A  A  A  9 

9  T  0 

A  A  A  4 

A  A  A  8 

»  A  A  3  rt 

RAO 

A  A  A  2 

F  9  B 

A  A  A  1  0 

HM  1 

A  A  A  1  1 

9  T  0 

A  A  A  2 

RAO 
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28 

18  0  2 

21 

18  5  6 

1  fl  0  9 

39 

18  0  9 

33 

1  Fl  1  3 

1  fl  H  9 

30 

1  8  H  9 

3  4 

1  f  5  6 

19  57 

31 

18  0  7 

33 

10  18 

19  4  5 

32 

184  5 

4  6 

1  9  4  B 

19  5  4 

33 

18  4  8 

60 

10  18 

1  A  7  4 

34 

187  4 

33 

19  7  7 

19  5  3 

35 

185  3 

4  6 

10  5  9 

I960 

36 

I860 

60 

18  18 

19  2  5 

3  7 

16  2  5 

20 

1  B  6  4 

19  5  7 

38 

13  5  5 

6  9 

18  6  2 

19  15 

39 

18  15 

34 

18  6  4 

1  B  6  •> 

4  0 

18  6  8 

33 

18  7? 

19  4  9 

4  1 

10  4  9 

46 

1  B  5  2 

1  »  1  7 

4  2 

1  8  S  2 

61 

19  2  4 

10  7  9 

4  3 

1  87  o 

34 

1  B  1  8 

19  5  7 

4  4 

1  8  S  7 

21 

18  5  6 

19  6  3 

4  5 

186  3 

24 

10  16 

19  19 

4  6 

1819 

39 

10  5  6 

19  6  5 

4  7 

18  6  5 

21 

19  2  0 

19  7  5 

4  3 

1875 

60 

1  B  1  6 

19  2  6 

4  9 

I  B  2  6 

32 

18  6  4 

19  4  1 

50 

184  1 

31 

18  6  4 

18  6  9 

51 

18  6  9 

60 

1  fl  3  0 

19  3  5 

5  2 

18  3  5 

32 

18  12 

19  4  0 

5  3 

16  4  0 

21 

18  3  0 

19  3  4 

54 

18  3  4 

61 

1  P  5  6 

19  6  6 

55 

16  6  6 

39 

18  2  0 

19  7  0 

56 

1  B  7  0 

31 

18  5  6 

19  7  6 

57 

18  7  6 

34 

1  B  3  0 

1  B  B  0 

58 

10  8  0 

21 

15  16 

1  fl  2  P 

59 

18  2  8 

67 

18  6  4 

1  fl  7  B 

to 

18  7  8 

20 

19  9  4 

19  17 

61 

18  3  7 

67 

1  B  1  6 

1°  3  1 

6  2 

163  1 

60 

90  0  2 

1  9  9  n 

63 

189  0 

34 

18  9  4 

1  fl  B  S 

64 

1  8  B  S 

3  3 

19  3  8 

19  9  1 

65 

183  1 

4  6 

10  17 

19  7  5 

66 

18  5  9 

60 

I  J  1  B 

IP  1  7 

t  7 

18  17 

60 

19  3  6 

1  B  9  1 

68 

18  9  1 

39 

18  6  4 

1  B  0  3 

69 

16  2  4 

10 

00  0  0 

00  5  1 

70 

18  1? 

30 

00  0  0 

00  5  1 

71 

1  8  3  B 

10 

00  0  0 

00  4  3 

72 

1  H6  2 

15 

70  7  9 

6  3  5  1 

7  3 

180  8 

7  A 

5  3  9  8 

16  5  0 

74 

18  7  3 

10 

00  0  0 

0  0  6  0 

7  5 

18  7  7 

10 

00  0  0 

0  0  4  n 
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APPENDIX  E 


Description  and  Explanation  of  the  IBM-650 
Computer  Program  Used  to  Calculate  Surface 
Heat  Flow. 


The  computer  program  was  written  to  calculate  the  heat  flow  out 
of  the  fuel  and  into  the  moderator.   The  neat  flow  out  of  the  fuel 
surface  is  given  .by 


(q/A)f(t)  =  -kf  j][    JlL       Bi  CD4  c°s  P^  +  ^  sln  P^) 


1=1  2a 


2   2 

n  n  a   ,.  2   2 

t  ,n  n  a. 


oo  ,JL.<*  at2  /  P        B        ( 2~} 


+E    2L  ;2         (  E 

n=l,3,5,-        —  1=1        "  "   "■    +  B.' 

.'  nna  ,£T4  ^1 


i  4L  (E-l) 


4   4    2 

1    TT 
/ 

16L 


s 

+ 


E 


q        a   A,    cosh   kL 
oo  j 


2    2  2    2 

n  tt  a   ,    ,    wn  tt  a  2. 


+ 


m       4L2  J        4L2 

x  t  A7  AT 

*       q  A.    e  N  V-"     cosh   (kL)    sinh   (V  — L  L) 

E  ;o: 21  - ,      c Ya       7X7   ia    - « •** - ] 

._.      k.   (ax     -  X.)  u  /j.  ,    x 

j-1        f  j  cosh(y  — L  L   ) 


and  the  heat  flow  into  the  moderator  is  given  by 

P 


(q/A)m    (t)  =   p   km  {  I  jE     B      (Di    cos   p    t  +   Ki    sin  p   t; 


1=1  2a 


2   2 

n  n   a      t  2    2 


/T2  rn  n  qN 

oo                 nrr.  4L               ,      p  D        v. — ^ — )  s  ,       T 

r                 2L*    e  /      f  Bi        4L2  ,    V-  qoo  a  A.    cosh  *L 

L.                     ,2  \      L  442             ./_,  22  22 

,    o    e           L  x     ,_ .  n  n   q           Q    2  ._  ,n  tt    q.  n  rr   a, 

n=l,3,5r*— -  i-1          a"  +  P,  J-1   km( — )+   \   )( — ) 

riTra  16L4  i  m     4L2  j         4L2 
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7 


-£      Q  F  A1  (     Ya        Slnh     V 

J  in     j  cosh  v'  — L  L 

a 


'ii 


where     D,   = 


CO 


sh  *;L      slnh  #;L   --cos  i/L      sin  i/L 


4  O  O  O  "3 

cos  i;L     Cosh  l/L  +  sin  i/L   slnh  i/L 
cosh  «/L  slnh  */L  +  cos  «;L   sin  J/L 


E  = 


cos'  Y;L     cosh  i/L  +  sin  «/L   slnh  V/L 

/FT 


and  o:  -     »/r 
-1    '2a 


The  equivalence  between  elements  of  the  algebraic  equation  and  the 
symbolic  logic  of  the  computer  program  is  shown  in  Table  E-l. 


Table  E-l.   Definition  of  symbolic  terms  of  the 
IBM- 650  computer  program  for  cal- 
culating surface  heat  flow. 


>  -  M 

i         7   2a 


™vi      -      cosh   /.L      sinh    J  L   -    cos    V.L      sin    V.  L 

ZMX1     -     i i i i _ 

cos      V  L     cosh      fcL  +    sin     fc.L      sinh     &  L 
A1i  =      ^iAi    (ZMX1i   cos  Pi11  -   ZMX2i    sin  Pi11    } 


aod       -       Q       a  A.   e     j 

A3P       =         co  j _ 

k(aK2  -  \.) 
J 


_       N  CO! 


/X  . 


_osh  kL     sinh  y  — L  L  ,    , 

A3SUM     -     7_a a  -   k  sinh  kL 


cosh  s/  — L  L 
a 


,     ^     cosjf  L   sin  I .  L  +  cosh#  L   sinhV.L 

'1  "      ■ =  E, 

cos     >tL  cosh2    VtL  +   sin2yiL  sinh2  ^L 
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Table  E-l   cont. 

A3.      =      (A3P.)    (A3SUM  ) 

CSHLL,     =   cosh    /_7  .    . 

COSLL.      =     cos     /    1   , 


ARG1        =      nn/2L 
n 

2    2  2 

ARG2        =      n   n   a/4L 
n 


2   2  4  4  2 

A2ST1        =      B,(2-H-a)/(   SJLA-     +  p    2   ) 
1      4L  16L 


2   2  2   2 

A2ST1.     =     q       a  A,   cosh  vh    I    k  ("  "  a  +  \.    )    i?-^^  +  a  k     ) 
j  °°  J  4L2  J  4L2 

2   2 

n    TT     TT     . 

nn  /T2 

2L        e     4L 

A2DDT        =     — — -= 

n  2    } 


P  s 

A2SUM       =  V       A2ST1        +  V  A2ST1 

1=1  j=l 


TERM   =  (A2DDT  )  (A2SUM  ) 
n  n        n 


The  Input  data  consists  of  the  heat  generation  and  surface 
temperature  parameters,  appropriate  material  constants,  half-thickness 
of  the  region  and  numerical  constants.   Table  E-2  lists  input  data 
required  for  the  program, 
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Table  E-2.  Input  Data  Required  for  Use  of  the  IBM-650 
Computer  Program  Used  to  Calculate  Surface 
Heat  Flow  Rates. 


Symbol 


Explanation 


Storage  Location 


ZERO  0.00  0164 

ONE  1.00  0656 

TWO  2.00  0820 

Pi  3.1415,9  0779 

FIFTY  50.00  0461 

CRIT  0.0001  0388 

ALPHA  Thermal  Diffusivity  0278 

KAPPA  Reciprocal  of  Thermal  Neutron  Diffusion  Length  0436 

in  Fuel 

Q00  Normalization  Factor  for  Heat*  Generation  0324 

KAY  Thermal  Conductivity  0581 

L  Half-thickness  of  Region  0186 

INDXM  No.  of  Terms,  Surface  Temperature  Fit  (OOOOOOOOxx)   0076 

NOLAM  No.  of  Terms,  Heat  Generation  Fit  (OOOOOOOOxx)  0456 

AYEJj  Amplitude  Parameter,  Heat  Generation  Fit  (0200  + 

LAMDAj  Exponential  Parameter,  Heat  Generation  Fit  (0020  + 

AMMMj  Amplitude  Parameter,  Surface  Temperature  Fit  (0100  + 

BTAAj  Period  Parameter,  Surface  Temperature  Fit  (0120  + 


J) 
J) 
J> 
J) 


The  output  from  this  program  is  punched  out  on  one  card  having  an 
eight  word  capacity,  one  word  consisting  of  10  digits  and  a  sign.  The 
form  of  the  output  is  shown  in  Table  E-3. 


Table  E-3.   Output  form  for  IBM-650  Computer 
Program  Used  to  Calculate  Surface 
Heat  Flow  Rates. 


WORD  1  WORD  2    WORD  3 


(q/A) 


n?     oo        s 

l     A1i   I(Term)n   I 


A3 


WORD  4   WORD  5  WORD  6   WORD  7   WORD  8 


j 


i=l 


nF  1,3,5 


j=l 


._ 

< 

—  _l 

<_l 

H-i- 

W-j  • 

Ouif 

F°  - 

0- 

O   «i    t> 

o 

u.  ot  ot 

_) 

.         (v 

\i 

O 

(E 

SbJ 

<=>  _l 

2  _l 

T  UJ 

10 

"  o 

H 

Z 

SS 

o 

en  to 

<j 

z 
o 
o 

2    . 

id 

o 

■   o 

o 

'~* 

tr 

%  • 

UJ 

ot  in 
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< 
+ 

fn    m  << 

<->  111  J* 


_J 

_l 

UK 

Q: 

U.X 

t- 

* 

*3 

2  UJ 

K 

in 

7 
1 

< 

r     l- 

CJ 

OT 

z 

F 

<t  o 

O 

o 
o 

(9 

s. 

z 

1- 

a  x 

u. 

o 

K 

UJ  2 

t| 

0> 

UJ    „ 

in  — 

t- 

z 

•>  • 

o 

in  in 

o 



+  + 

\       a. 

o-h 

1         o 

r-h 

_i 

UJ  UJ 

Ei- 

u°"  ,7 

.m  L 

.        V    V 

£  ot  ot 

<M     ■    ■ 

«S  -"-" 

,hl- 

E  </)  in 

5  C\J<\J 

li  <  < 

.  2S       2 

"  £  tr  •-  rr 

2o£o£ 
<r  U-T  -i  . 

UJ    •+  o  + 

-S31 3 

o 

cUO  cO 

H 

lONini 

en 

c 

fo\_ 

«      IS 

o      m 

\zj~ 

2 

_J        UJ 

tr 

O       h- 

UJ 

I    c 

|  tr  tr 

E 

J= 

Q 

O  UJ  O 

O 

li 

a 

u.  t-  u. 

o 

< 

UJ 


a 

z 

UJ 

a. 
a. 

< 
i 

< 
rr 

< 
a 
o 
e> 
o 
_l 


■0»- 

"  2 

o  cB 

o 

OT12 

IE 

o   . 

< 

<?  e 

h 

cn 

/o^ 

ii 

W 

-  tB 

<\4  — 

E 

X  ~T 

o 

2<0. 

N27 

o.i' 

o  c 

—  in  » 

a.  ~ 

m 

Eg 

z 

■o 

a. 

£  in 

ir> 

a. 

o  o- 

UJ 

tr  o 

o 

u.  o< 

aj 

_i 

Form  COSLL 
Form    A3SUMj 
Form    A3| 
Set  i'j-l 

o 

UJ 

z 

2 
_l  3 
JOT 

inio 
w«<7 

°         ir 

6  !  It 

OT 

o 

Cl 

+  °    n    (ttt 
Z  Z  f  <  < 

z  _<  E  E 


o 

UJ        t£ 
Z       UJ 

■  Z  a 

<°Q 

I'o 
5  zi 
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BLR 

0  20  0 

0260 

Rl  R 

18  0  0 

189  9 

8  r  n 

A  R  T  A  H 

18  00 

«YN 

START 

199  9 

8  Y  N 

A  Y  E  J 

0  3  00 

8  Y  N 

L  A  M  0  A 

0  330 

8  Y  N 

A  3 

0  34  0 

BLR 

19  5  1 

1960 

BLR 

19  7  7 

1  98  4 

BLR 

0  mo 

0160 

8  Y  N 

A  U  U  U 

0  100 

»  Y  N 

B  T  A  A 

0  130 

8  Y  N 

ALPHA 

0  3  7  8 

8  Y  N 

KAPPA 

0  4  36 

3  Y  N 

L 

0  18  6 

8  Y  N 

000 

0  3  3  4 

3  Y  H 

KAY 

0  58  1 

8  Y  N 

1  H  [1  I  » 

00  7  6 

8  Y  N 

N  0  L  A  M 

0  4  56 

8  Y  N 

A  1 

0140 

t  00£  • 

3  TO 

A  A  A  1 

8  I  U 

A  A  A  1  4 

R  A  U 

A  A  A  1  6 

F  A  M 

A  A  A  1  4 

S  T  U 

A  A  A  2 

LOO 

A  A  A  3 

8  T  0 

A  A  A  4 

A  A  A  8 

lllfl 

R  A  U 

A  A  A  2 

FSB 

A  A  A  5 

BU  1 

A  A  A  6 

S  T  U 

A  A  A  2 

R  A  U 

A  A  A  4 

F  M  P 

A  A  A  7 

3  T  U 

A  A  A  4 

A  A  A  8 

A  A  A  6 

R  A  U 

A  A  A  2 

FSB 

A  A  A  3 

BM  1 

A  A  A  2  8 

8  TU 

A  A  A  2 

RAO 

A  A  A  4 

FUP 

A  A  A  9 

S  T  U 

A  A  A  4 

A  A  A  6 

A  *  A  38 

R  A  U 

A  A  A  2 

F  S  R 

A  A  A  1  0 

BM  1 

A  A  A  1  1 

8  T  0 

A  A  A  2 

R  A  U 

•  A  A  4 

F  M  P 

A  A  A  1  2 

S  T  0 

A  A  A  4 

A  A  A  38 

A  A  Al  1 

RAO 

A  A  A  2 

LOO 

A  A  A  1  7 

FUP 

A  A  A  4 

S  TO 

A  A  A  1  3 

R  A  U 

A  A  A  1  4 

8  M  1 

A  A  A  1  5 

RAO 

A  A  A  1  3 

A  A  A  1 

A  A  A15 

RAO 

A  A  A  3 

F  0  V 

A  A  A  1  3 

A  A  A  1 

A  A  A  17 

8  T  0 

A  A  A  1  8 

RAO 

A  A  A  3 

FAD 

A  A  A  3 

3  T  0 

A  A  A  1  9 

LOO 

A  A  A  27 

.  8  T  0 

A  A  A  20 

8  T  0 

A  A  A  2  1 

'RAO 

A  A  A  2 

FUP 

A  A  A  2 

8  T  0 

A  A  A  2  3 

A  A  A  23 

A  A  A2  2 

F  D  V 

A  A  A  2  1 

S  T  0 

A  A  A  24 

F  A  0 

A  A  A  1  9 

S  TU 

A  A  A  1  9 

RAO 

A  A  A  24 

F  D  V 

A  A  A  1  9 

FSB 

A  A  A  2  5 

BU  1 

A  A  A  2  6 

RAO 

A  A  A  1  9 

A  A  A  1  8 

A  A  A36 

RAO 

A  A  A  20 

F  A  0 

A  A  A  3 

8  T  U 

A  A  A  20 

FU  P 

A  A  A  2  1 

8  T  0 

A  A  A  21 

RAO 

A  A  A  2  3 

FUP 

A  A  A  2 

8  TU 

A  A  A  23 

A  A  A  32 

A  A  A  3 

10 

00  0  0 

00  5  1 

A  A  A  5 

5  0 

00  0  0 

00  5  1 

A  A  A  7 

14 

8  4  10 

00  5  3 

A  A  A  9 

37 

18  30 

00  5  1 

A  A  Al  0 

20 

00  00 

00  5  0 

A  A  A  1  2 

13 

2  14  0 

00  5  1 

A  A  A  1  6 

00 

OO  00 

00  00 

A  A  A2  5 

70 

0  000 

0047 

A  A  A  3  7 

30 

0  0  n  0 

00  5  1 

I   U  X  3  B 

8  T  0 

A  B  C  0  1 

RAO 

A  R  G  GL 

LOO 

C09HX 

8  T  0 

C  H  A  AL 

RAO 

A  R  C  G  L 

LOO 

E  0  0  C  Ft 

8  T  0 

C  8  A  AL 

RAO 

A  R  G  G  L 

LOO 

8  1  N  H  X 

8  T  0 

8H  A  AL 

R  A  U 

A  R  G  GL 

LOO 

EOOSH 

8  T  U 

3  N  A  AL 

RAO 

C  H  A  AL 

FUP 

C  3  A  AL 

S  T  U 

F  CL 

F  M  P 

C  H  A  AL 

FUP 

C  8  A  AL 

8  T  0 

0  E  1  1  H 

R  A  U 

8  N  A  A  L 

FUP 

8  H  A  A  L 

8  T  II 

F  3L 

1 

00  00 

00 

00  0  0 

00  0  0 

2 

0  0  0  0 

0  0 

00  0  0 

0  0  0  0 

3 

0  0  0  0 

00 

00  0  0 

0  00  0 

4 

0  0  0  0 

0  0 

00  0  0 

0  00  0 

5 

0  0  0  0 

00 

00  0  0 

00  0  0 

6 

00  0  0 

00 

00  0  0 

0  00  0 

7 

00  0  0 

00 

00  0  0 

000  0 

8 

00  0  0 

00 

00  0  0 

000  0 

9 

0  0  0  0 

00 

00  0  0 

00  0  0 

1  0 

00  0  0 

00 

00  0  0 

0  0  0  0 

1  1 

0  0  0  0 

00 

00  0  0 

0  0  0  0 

13 

0  0  0  0 

00 

00  0  0 

J  0  0  0 

1  3 

0000 

00 

00  0  0 

00  0  0 

1  4 

0  0  0  0 

00 

00  0  0 

00  0  0 

1  5 

00  00 

00 

00  0  0 

0  0  0  0 

1  6 

00  0  0 

00 

00  0  0 

OOO  0 

1  7 

00  0  0 

00 

00  0  0 

0  00  0 

1  8 

0  0  0  0 

00 

00  0  0 

00  0  0 

1  9 

00  0  0 

00 

00  0  0 

0  0  0  0 

SO 

0  0  0  0 

00 

00  0  0 

00  0  0 

21 

0000 

2  4 

00  0  3 

00  0  6 

32 

P  0  0  6 

31 

00  10 

00  13 

2  3 

00  13 

60 

00  16 

00  2  1 

2  4 

0  0  2  1 

37 

00  1  0 

00  3  7 

2  5 

0  0  3  7 

21 

00  4  2 

0  0  4  5 

26 

0  0  4  5 

6  9 

00  4  8 

0  00  1 

27 

00  0  1 

2  4 

00  0  4 

00  0  7 

28 

0  0  0  7 

60 

00  4  2 

00  4  7 

2  9 

00  4  7 

3  3 

00  5  0 

00  2  7 

30 

0  0  2  7 

4  6 

00  3  0 

00  3  1 

31 

00  3  1 

21 

00  4  2 

00  9  5 

33 

0  0  9  5 

60 

00  0  4 

0  0  0  9 

3  3 

0  0  0  9 

3  9 

00  1  2 

0  0  6  2 
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The  safe  operation  of  nuclear  reactors  i"s  imperative  if  there  is  to 
be  increased  engineering  application  of  these  systems.   Transient  reactor 
experiments,  such  as  the  SPERT  tests,  have  demonstrated  that  thermal, 
heterogeneous,  liquid-moderated  reactor  systems  will  safely  shut  them- 
selves down  following  step  and  ramp  insertions  of  limited  amounts  of 
excess  reactivity.   It  is  important  that  a  model  based  on  the  nuclear, 
thermodynamic  and  hydrodynamic  properties  of  the  reactor  system  be 
developed  to  explain  this  phenomena  so  that  it  can  be  used  in  the  design 
of  new  systems. 

Equations  for  the  fine  structure  of  the  temperature  distribution  in 
a  unit  cell  of  a  heterogeneous  reactor  during  a  transient  burst  were 
derived  based  on  the  known  power  and  fuel  surface  temperature  distributions. 
A  model  based  on  recognized  shutdown  effects  was  developed  to  calculate 
the  excess  reactivity  during  a  transient  using  the  temperature  distribu- 
tions to  define  the  deposition  of  energy.   The  calculated  excess  reactivities 
show  this  model  to  be  satisfactory.   The  effect  on  reactivity  due  to  steam 
formation  required  one  empirical  parameter  which  can  probably  be  removed 
when  a  greater  knowledge  of  transient  boiling  is  available. 


Date  Due 

• 

